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PREFACE 


At its thirteenth session, held in Paris in October-November, 1964, 
the General Conference of the United Nations Educational, Scientific and 
Cultural Organization (UNESCO) set up the Co-ordinating Council of the 
International Hydrological Decade. The work of the Council lead to an 
international program in the field of hydrology. The overall objectives 
were to accelerate the study of water resources and the regimen of 
waters with a view to their rational management in the interest of 
mankind, to make known the needs for hydrological research and education 
in all countries, and to improve each country's ability to evaluate its 
water resources and use them to the best advantage. The program thus 
focused on science, but gave strong consideration to utilitarian factors. 

As part of the contribution to the International Hydrological 
Decade program, the Province of Ontario investigated the water resources 
and physical conditions in five drainage basins in Southern Ontario, 
each being representative of different conditions common in the province. 
The hilly moraine and mild karst topography of the Blue Springs Creek 
basin was selected as one of the representative areas. This report 
describes the hydrogeologic characteristics of the basin as deduced from 
data collected during the International Hydrological Decade, 1965-1974. 
The study was completed by the Water Resources Branch of the Ontario 


Ministry of the Environment. 
: lay be Jer 
at 


G.H. Mills, Director 
Water Resources Branch 
Toronto, February 1978 
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ABSTRACT 


A study of the occurrence, quantity and quality of the ground water 
was undertaken in the Blue Springs Creek basin by the Ontario Ministry 
of the Environment, as part of the International Hydrological Decade 
(IHD) Representative Basin program. The Blue Springs Creek basin was 
selected for study to represent a physiographic region in Southern 
- Ontario, characterized by shallow overburden and karst limestone and 
dolomite bedrock topography. These conditions are generally found along 
the top of the Niagara Escarpment. 

Field work for this study included the investigation of the sur- 
ficial geology, the installation of 16 observation wells to aid in the 
investigation of the subsurface geology and ground water movement, 
streamflow monitoring, and chemical sampling of surface water and ground 
water systems. 

The major aquifer found in the basin is the Amabel Dolomite Formation; 
it averages 60 feet (18 metres) in thickness, reaches a maximum thickness 
of 100 feet (30 metres) and outcrops over approximately 10 per cent of 
the basin area. The overburden material, covering the remaining 90 per 
cent, forms a minor aquifer, composed largely of till and kame moraine 
deposits with a thickness of almost 150 feet (45 metres) in the northern 
portion of the basin. 

The Amabel Formation in the basin was found to have a median 
transmissivity of 3,200 IGPD/ft (550 mm /s;) similar to that of the 
Overburden aquifer at 3,500 IGPD/ft (600 mm/s) . The transmissivity 
values found for shallow wells in the Amabel Formation were generally 
higher than those found in deeper wells in the same formation. 

Baseflow in the basin was found to average approximately 68 per 
cent of the total streamflow during the water years of 1966-69, and was 
equivalent to an annual ground water runoff rate of 9.8 inches/yr (250 
mm/yr). The average specific yield was found to vary from 7.2 to 12.7 
per cent in various sub-basins of the Blue Springs Creek basin. 

Two finite difference ground water models were used to simulate the 
steady state and dynamic flow patterns in the basin, and to improve or 
verify the estimates of transmissivity, coefficient of storage, recharge 
and baseflow obtained by conventional techniques. Good representation 
of the piezometric surface was obtained with the model if the regional 
bedrock transmissivities used in the model were 1.5 times the median 
transmissivities found at wells, the storage coefficients ranged from 4 
to 16 per cent and the recharge rates varied from 7 to 15 inches per 
year (180 to 380 mm per year) in different areas of the basin. 

In examining the hydrochemistry of the flow systems, the majority 
of the ground water samples taken from bedrock wells indicated that the 
water was of the calcium-magnesium and bicarbonate-chloride sulphate 
type, saturated with respect to calcite and undersaturated with respect 
to dolomite. 

The calcium to magnesium ratios of the water samples ranged from 
1.31 to 4.61 with a mean of 2.10 in comparison with a mean ratio of 1.09 
in the bedrock material and 2.15 in the overburden material. The higher 
ratios in the ground water tended to be in the recharge areas, and the 
lower ratios in the discharge areas. Using these results the generalized 
ground water flow patterns were determined, and broadly confirmed the 
piezometric flow patterns. 


INTRODUCTION 
SCOPE OF THE INVESTIGATION 


The study of the hydrogeology of the Blue Springs Creek basin is 
one of several detailed water resources investigations being completed 
by the Ontario Ministry of the Environment, as part of the International 
Hydrological Decade (IHD) Representative Basin program (Figure 1). The 
University of Guelph studied the surface water hydrology as an IHD co- 
operative project with’the Ministry of the Environment. 

The major areas of study described in this report pertain to investi- 
gations of the geology, hydrogeology and hydrochemistry of the Blue 
Springs Creek basin. Computer models were applied to simulate the 
hydrogeologic conditions. 

The study included field investigations of the geology, hydrogeology 
and hydrochemistry. Water well records maintained by the Ministry of 
the Environment were used to characterize the aquifers, water levels and 
subsurface geology in the region. The streamflows in the basin were 
monitored by three stage recorders maintained by the Ministry of the 
Environment and Environment Canada, as well as five staff gauges which 
were read about once a month. Ground water levels were monitored by four 
recorders set up on observation wells, and by approximately two monthly 
readings taken at four other wells. Aerial photographs obtained from 
the Ministry of Natural Resources were used to supplement field studies. 


OBJECTIVES 


The basic objectives of the study are summarized below. They are 
based on the 'Guidelines for Research Basin Studies' (1966), as compiled 
by the Canadian National Committee for the International Hydrological 
Decade: 

iA To study the surficial and subsurface geology of the basin and 

its surrounding area. 

24 To study the geologic conditions under which various flow 

systems operate in the basin. 

aa To determine the hydrogeologic parameters of the aquifers in 

the basin. 

4. To identify ground water flow systems in the basin. 

55 To delineate principal areas of recharge and discharge in the 


basin. 

6's To evaluate the total ground water discharge. 

Ta To determine the basin yield. 

8. To estimate the regional hydrogeologic parameters of the 
bedrock aquifer using a ground water model. 

9. Ta determine the relationship between water quality and water 
movement. 


DESCRIPTION OF THE BASIN 


The Blue Springs Creek basin comprises an area of approximately 30 
Square miles (78 km*), bounded by latitudes 43° 30' and 43° 4E° Nand 
longitudes g0° 00' and 80° 10' W. It contains part of Wellington and 
Halton counties. : 

The basin is almost entirely rural, including only a few small 
settlements such as the villages of Knatchbull and Crewsons Corners; the 
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larger towns of Acton, Eden Mills and Rockwood lie just outside the 
drainage boundary (Map 1). When considering the overall hydrologic 
regime of the basin, it can be considered as relatively unaffected by 
human activity. 

The topography of the Blue Springs Creek basin (as shown on Map 
is characterized by broad rolling hills, interrupted by small and 
prominent ridges and deeply incised stream valleys. The maximum elev: 
tion in the basin is about 1400 feet (425m) above mean sea level, in tt 
northern portion of the basin. The lowest elevation is about 1040 feet 
(315m), at the downstream end of the basin, giving a total relief of 
approximately 350 feet (110m). 

The basin is dissected by the perennial Blue Springs Creek, which 
is a tributary of the Eramosa River in the Grand River watershed. Its 
length is about eight miles (13 km), with an average gradient of approx 
mately 45 feet per mile (8m per km). The creek has its main source in 
the morainic hills near Acton and is also supplied by a number of spring 
in the basin (Dickinson and Whiteley, 1970). The creek derives its name 
from a large spring in the centre of the basin. 

The basin contains several intermittent streams (Map 1), some of 
which do not reach the main channel but sink into local depressions. 
These springs and sinking streams are further discussed in Section 
(2, 224). 

The climatic and hydrologic data collected in the period 1967-1971 
were used for much of the hydrologic analysis, as they formed the most 
complete set of data that was available when the analysis was carried 
out. These five years were slightly wetter than the historic period as 
recorded at nearby long term stations. Thus the mean precipitation at 
Guelph during 1967-1971 was 33.7 inches compared to the 30 year (1941- 
1970) mean of 32.7 inches. The mean streamflow on the Speed River below 
Guelph during 1967-1971 was 212 cfs compared to 197 cfs for the 24 years 
of record (1950-1973). As the mean conditions during the years 1967- 
1971 were only slightly wetter than the long term mean values (3% and 8% 
higher for the precipitation and streamflow, respectively), the hydro- 
logical conditions during this five year period were considered to be 
fairly typical of the long term conditions. 


PREVIOUS INVESTIGATIONS 


The Paleozoic rocks underlying the area were previously mapped by 
Williams (1919), Caley (1941), Bolton (1957), Sanford and Quilliam 
(1959), Sanford (1961 and 1969), Liberty (1969), Hewitt (1972), and 
Telford (1976). Hydrogeologic information on the Paleozoic rocks under- 
lying the area is contained in water well records maintained by the 
Ministry of the Environment, and oil and gas well records maintained by 
the Petroleum Resources Section, Ministry of Natural Resources. Hewitt 
(1964 a, b) has also discussed the limestones of the area for use as 
building materials and as industrial minerals. 

The Pleistocene features and deposits of the area have been mapped 
by Chapman and Putnam (1966) and by Karrow (1968). 

Agricultural soils have been mapped by Hoffman and Mathews (1963) 
and by Gillespie et al (1971). 

The study of watershed areas contributing to runoff in the Blue 
Springs Creek basin, by Dickinson and Whiteley (1970), also provided 
significant background data for the study. 
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GEOLOGY 


The characteristic features of the bedrock and surficial geology of 
the Blue Springs Creek basin are discussed in this chapter. A summary 
of the various units and their position in the stratigraphic column is 
given in Table 1. The distribution of surficial geologic units is shown 
on Map 2, and cross-sections of the basin geology are illustrated in 
Figure 3 and Figure 5. 


BEDROCK GEOLOGY 
Regional 


The consolidated bedrock formations in southwestern Ontario consist 
of Paleozoic sedimentary rocks of Ordovician, Silurian and Devonian Age, 
resting on the Precambrian basement (Figure 2 in pocket). The Paleozoic 
rocks dip regionally to the southwest at about 25 feet to the mile; the 
Oldest rocks of Precambrian Age outcrop to the northeast and the younger 
rocks of Devonian Age outcrop in the southwestern part of Ontario. In 
the area around Blue Springs Creek, the outcropping rocks are dolomites 
of the Amabel Formation of Middle Silurian Age (Figure 2); older rocks 
of Lower Silurian and Ordovician Age have been found in the subsurface 
in a number of oil and gas wells drilled in the area (Figure 3 in pocket). 

The Silurian rocks in Southern Ontario undergo a facies change 
north of Hamilton, which has been attributed either to a structural 
barrier called the 'Algonquin Arch' (Roliff, 1954), or to the presence 
of two major depositional areas during Silurian times in Southern Ontario 
(Bolton, 1957). The Silurian rocks contain many carbonate beds, which 
to the north, in the Bruce Peninsula area, (Figure 2) are more dolomitic, 
while to the south, in the Niagara Peninsula, they are more calcareous. 
In the Blue Springs Creek area, the Clinton Group (Table 1) has been 
correlated with the more calcareous, southerly facies of the Niagara 
Peninsula, while the overlying Albemarle Group is more closely cor- 
related with the dolomitic northerly facies (Bolton, 1957). 

In the Blue Springs Creek area, minor warping of the Ordovician and 
Silurian beds occurs (Liberty, 1969). Weber (1960) considers that the 
Eramosa River near Rockwood, (Map 1), runs along an anticlinal axis in 
the Amabel Formation. The southern part of Blue Springs Creek runs 
along a synclinal axis in the same formation (Liberty, 1969), but this 
fold does not extend into the northern portion of the basin, where the 
Amabel is relatively unfolded. 


Blue Springs Creek Basin 


In the Blue Springs Creek basin, the outcropping bedrock beds are 
members of the Middle Silurian Amabel Formation (see Table 1 and Map 2). 
The Upper Ordovician and Lower Silurian beds, however, also play a 
significant part in the hydrogeologic characteristics of the basin . 
(Figure 5 in pocket). In addition, these beds form the source material 
from which much of the overlying glacial deposits were derived. 

The Ordovician beds, resting unconformably on the Precambrian . 
basement, have a total thickness of about 1900 feet in the Blue Springs 
Creek area (Table 1). The lower 1450 feet of the Ordovician beds are 


composed of various sedimentary rocks, and are covered by the Upper 
Ordovician Queenston shale. This is a soft red impervious shale about 
460 feet thick. The Ordovician beds outcrop to the east of the basin 
(Figure 2) and therefore form the source material for much of the 
glacial deposits in the basin, as the Wisconsinan ice mass that gla- 
ciated the basin moved from the southeast (see Section 2.2.1). 

The Queenston shale is covered by the Silurian Cataract and Clinton 
Groups, which together consist of about 85 feet of shales and dolomites 
(Table 1). None of these beds outcrop in the basin, but they have been 
observed during the drilling of wells in the area (Figure 3 in pocket). 
Overlying the Clinton Group is the Middle Silurian Albermarle Group, 
comprising the Amabel and Guelph Formations. The Amabel Formation has 
been divided into four members, three of which are present in the Blue 
Springs Creek basin (Bolton, 1957, Sanford, 1969). The Lions Head 
Member, which is the lowest member of the Amabel group, occurs in the 
Bruce Peninsula but is not present in the Blue Springs Creek basin. 

The lowest member of the Amabel Formation present in the Blue 
Springs Creek basin is the Colpoy Bay Member (see Table 1), which is a 
buff weathered, massive, porous, fine grained to dense, white to blue 
grey dolomite (Bolton, 1957), lying unconformably on the Clinton Group 
rocks. The overlying Wiarton Member is a similar dolomite but is light 
grey with blue grey mottling. Sanford (1969) does not differentiate 
between the Colpoy Bay and Wiarton members but groups them together as 
the Wiarton Formation. Bolton (1957), however, places the Colpoy Bay- 
Wiarton contact at the base of the mottling which is discernible at 
various exposures in the general area. Both the Wiarton and Colpoy Bay 
members contain numerous bioherm reefs, which are usually more porous 
than the surrounding rocks, and are important cave-forming structures. 
The Wiarton Member was identified in outcrops along the channel of Blue 
Springs Creek (Map 2). The Colpoy Bay Member, however, was not identi- 
fied in the basin, but it outcrops or is present just below the stream- 
bed along the channel of Blue Springs Creek, as indicated in Figure 3. 
The overlying Eramosa Member is a bituminous, thin bedded, dark brown, 
fine grained to sugary dolomite (Bolton, 1957), which outcrops in the 
extreme western part of the basin (Sanford, 1969). The overlying Guelph 
dolomite, which is also in the Albermarle Group, is not present in the 
basin, but outcrops to the west of the area. 


Bedrock Topography 


Elevation contours of the bedrock surface have been drawn utilizing 
data from records of water wells drilled in the study area. The locations 
of these wells are shown on Map 3 and the bedrock surface elevations are 
shown on Map 4. 

The bedrock surface in the Blue Springs Creek basin has a maximum 
relief of almost 200 feet (60m), which appears more subdued than that of 
the present land surface. The lowest bedrock elevations are found‘in 
the vicinity of Eden Mills where elevations are less than 1050 feet 
(320m) above sea level; the highest elevations, in excess of 1225 feet 
(370m) above sea level, occur in the northern portion of the basin. 

The most pronounced feature of the bedrock topography is the 
bedrock valley which is now occupied by the main channel of Blue Springs 
Creek. This channel was formerly a glacial spillway in front of the 
Galt and Moffat moraines (Section 2.2.1); further deepening has sub- 


sequently taken place through recent erosion of the Blue Springs Creek 
channel. 


TABLE 1. STRATIGRAPHIC COLUMN OF GEOLOGICAL UNITS IN THE BLUE SPRINGS 
CREEK AREA. AFTER BOLTON (1957) AND SANFORD (1969) 


EPOCH GROUP FORMATION MEMBER CHARACTER OF APPROXIMATE 
MATERIAL THICKNESS 
IN BASIN 
(Feet) (metres) 
Recent Alluvium gravel,sand,silt variable 
Swamp or muck,marl,peat variable 
Marsh 
Pleis- Late Wis- silty sand till, 0-150 (0-50) 
tocene consinan kame,esker & out- 
Glacial wash sand & gravel 
Driot . 
Middle Alber- Guelph (not present in the Blue Springs O 
Sil- marle Creek basin) 
urian et 
Amabel Eramosa bituminous brown 0-30 (0-9) 
dolomite ‘2 
Wiarton massive, light grey 30 (10) 
porous, mottled, 
dolomite re 
Colpoy massive, white to SO (15) 
Bay blue grey porous, 
dolomite ah 
Lions (not present in the Blue @) 


Head Springs Creek basin) 


Clinton Reynales massive light grey He. C5)) 
dolomite 
Lower | Cataract Cabot grey shale 40* (12) 
Sil- Head Er 
urian 
Manitou- argillaceous 20* (6) 
En dolomitic 
limestone ae 
Whirlpool massive sandstone HO Seco) - 
Upper Queenston red shale 460* (140) 
Ordo- 
vician = 
Ordo- Various various 1450* (450) 


vician 
e * 2 
Precam- Various various E 


brian 


*thickness from gas well, Nassagaweya Fh peas 
Cone) sl LL. 


Karst Features 


Karst landforms consist of features which have been developed 
predominantly by solution. Typical karst features are developed in 
limestones, but can also occur in gypsum, rocksalt, dolomite and other 
water soluble rocks. Dolomite rocks are less soluble in water than 
limestones (Sweeting, 1973) and therefore have less rapid karst develop- 
ment, as found in the Amabel dolomites in the Blue Springs Creek basin. 

One feature of karst areas is the extensive development of under- 
ground drainage, which in mature karst leads to the complete lack of 
perennial surface streams. The initiation of such underground drainage 
occurs in bedding planes and fractures, which become enlarged due to 
progressive solution of the rock. The initial development of the 
solution channels is almost always under phreatic conditions (below the 


water surface); however, as the openings enlarge, development can 
continue under phreatic or vadose (free air to water surface) conditions 
(Ford, 1972). As these conduits enlarge sufficiently, they can even- 


tually form cave structures. 

In the Blue Springs Creek basin, no extensive caves have been 
found; however, small caves do occur nearby at Rockwood, (Photograph 1) 
one mile to the west of the basin (Map 1). These have been described by 
Paton (1889), Weber (1960), Kershaw (1973) and Cowell and Gregor (1973). 

A number of theories for the origin of the Rockwood caves have been 
proposed. Weber (1960) thought their formation was the result of an . 
anticline at Rockwood, although this theory has not been supported by 
other writers. In many other karst areas the solution channels form in 
the more calcitic and purer members of a limestone/dolomite sequence 
(Rauch and White, 1970, Sweeting, 1973). Bradley (1968) considered that 
the Rockwood caves were formed in the bioherms of the Amabel Formation, 
which he considered were more calcitic. Cowell and Gregor (1973), 
however, have analyzed the composition of the dolomitic bedrock and the 
bioherms in the Rockwood gorge for calcite and dolimite by atomic 
absorption spectrophotometry. They found no significant difference 
between the compositions of the bioherms and the adjacent bedded dolo- 
mites (the mean molar calcium to magnesium ratio of two samples of the 
bioherms was found to be 1.02, and in three samples of bedded dolomite 
was 1.05. Ideally, pure dolomite (CaMg (CO3)5) would have a molar 


calcium to magnesium ratio of 1.0). The composition of the bioherms 
therefore, appears to have played no significant part in the cave 
development. 

The Rockwood caves occur near the bioherms of the Amabel Formation 
(Cowell and Gregor, 1973); however, no caves were found in the bioherm 
structures themselves, as stated by Karrow (1968) and Weber (1960). The 


cave structures appear to have developed more readily in the dolomite 
beds near or at the base of the bioherms, where structuring has dis- 
turbed the bedding and increased the porosity of the dolomite. 

Cave orientation in the dolomite formation at Rockwood also has 
been discussed by Weber (1960) and Cowell and Gregor (1973). The latter 
authors suggest that the orientation is largely controlled by the 
bioherms. The predominant reef direction (325° from north) (Karrow, 
1968) does play a significant part in the orientation of the cave 
development. This control would also have an influence on the permea- 
bility of the dolomite, as conduits would tend to be more open and 
permeable in the direction of the bioherm reefs. 

Other major features that characterize typical karst areas are 


Photograph 1: Caves and reSurgences found at the conservation 
area in the vicinity of Rockwood. 


Photograph 2: Sinkholes are found at ground surface where 
dolomite bedrock is exposed or covered by 
only a thin layer of till. 


sinkholes and springs. Both types of features can be observed to a 
limited extent in the Blue Springs Creek basin. 

Sinkholes appear as actual openings in the ground surface, usually 
connected to vertical solution channels that lead to other underground 
drainage systems. The sinkholes observed in the study area appear 
relatively small in size, occurring locally in areas where the dolomite 
bedrock is exposed or covered by only a thin layer of till (Photograph 
2). The 'sinking' phenomena can also be observed along various tribu- 
tary stream channels where streams actually flow into local depressions 
and disappear through the soil and overburden materials. This usually 
occurs in areas where the dolomite bedrock is close to the surface. 
Surface channels downstream of these sinks are usually very poorly 
defined as flows very seldom reach the main stream. Several such 
streams have been identified in the basin (Map 1). Some are inter- 
mittent streams, only flowing during wet periods, while others appear 
perennial in their course. 

Two major types of karst springs or resurgences can be found in the 
area, those where water discharges from vadose channels and the vauclu- 
Sian type, where the water issues from depth out of phreatic channels 
(Sweeting, 1973). Photograph 1 shows some springs at Rockwood of the 
vadose or free flowing type, while the Blue Spring itself, the location 
of which is shown on Map 1, is of the vauclusian type, the water rising 
through a mantle of till. Hydrologically, the two types of springs are 
significantly different (White, 1969). The free flowing springs, as at 
Rockwood, act like surface streams, exhibiting pronounced floods and 
droughts; whereas the vauclusian springs are fed from a relatively large 
ground water reservoir, tending to smooth out any peaky flows. In the 
Blue Springs Creek basin, the springs shown on Map 1 are largely fed 
from phreatic channels and so would tend to have less peaky flows than 
the surface streams. 

In typical karst areas the aquifer characteristics, such as trans- 
missivities and storage coefficients, are extremely heterogeneous, with 
the major portion of the ground water flow taking place along solution 


channels or caves. In mature karst, the flow of ground water can be 
represented by flow through a pipe network and not through a granular 
medium, as is usual for clastic aquifers (White and Longyear, 1962). In 


attempting to develop piezometric maps for a mature karst area, extreme 
ground water gradients would likely be found, as adjacent wells may tap 
water in channels which could have piezometric levels many feet apart. 
Such a map would therefore contain many discontinuities and many spurious 
points where the static level would be anomalously low or high. MIn- 
spection of the piezometric map (Map 7) of the Blue Springs Creek basin, 
shows that the piezometric surface is well behaved and that there are 
very few spurious points. These findings tend to indicate that the 
karst in Blue Springs Creek is not well developed and that the flow in 
the dolomite must be taking place through the joints and bedding planes. 
The ground water flow pattern, therefore, has not been significantly 
affected by the karst development. 

In mature karst areas, the drainage basin (surface and subsurface) 
is often considerably different in size from the topographic basin, due 
to water crossing topographic divides in underground solutional channels 
(White and Schmidt, 1966). In the Blue Springs Creek basin, some small 
differences have been observed between the topographic and the drainage 
divides (Map 7), where ground water moves across the topographic divide 
in the dolomitic joints and bedding planes. 
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SURFICIAL GEOLOGY 


The overburden deposits in the Blue Springs Creek basin consist of 
glacial drift of Pleistocene age with minor amounts of alluvial and 
swamp deposits of Recent age. The Pleistocene deposits consist mainly 
(Map 2). The thickness of the overburden material (Map 5) ranges from a 
few feet in the main creek valley, gradually increasing towards the 
basin divide, to a maximum thickness of almost 150 feet in the northern 
part of the basin (see also Figure 5). 


Glacial History 


The glacial features that are present in the Blue Springs Creek 
basin are the result of the late Wisconsin ice advance, which retreated 
about 12,000 years ago (Karrow, 1963). The basin was overrun by several 
pre-Wisconsin ice movements; however, the Wisconsin ice advance has 
obliterated evidence of these former glaciations in the basin. 

During the Wisconsin glaciation, the ice flowed from the present 
Lake Ontario basin and crossed the Blue Springs Creek area in a west- 
northwest direction, as is evident from the orientation of drumlins 
(Chapman and Putnam, 1966). The ice laid down till plains of the 
Wentworth till, which was derived from the Silurian and Ordovician 
formations, outcrops of which can presently be found to the east of the 
area. Drumlins were formed under the ice, and occur as isolated hills 
on the till plains. As the ice retreated, pauses or possible slight 
readvances deposited the Paris, and later the Galt and Moffat moraines 
(Karrow, 1968). Meltwaters flowed down to the southwest, in front of 
the Paris moraine, and carved out a prominant spillway, part of which is 
now occupied by the Eramosa River. Meltwaters in front of the Galt 
moraine flowed southwest from north of Acton and then broke through the 
Paris moraine at Eden Mills; meltwaters in front of the Moffat moraine 
flowed across the Galt moraine about two miles to the southwest of Acton 
(Map 1). The main channel of Blue Springs Creek now occupies these 
spillways. 

As the ice retreated, kames, eskers and outwash plains and terraces 
were deposited by the sediment laden waters flowing off the ice masses. 
All these deposits are present in the basin (Map 2) and are discussed in 
BECTION 2.2. 3. 


Wentworth Till 


Glacial till can be described as a nonstratified sediment, carried 
and deposited by glacier ice. Grain sizes in tills vary from clay to 
boulder size and can be mixed in various proportions. The Blue Springs 
Creek basin contains the sandy Wentworth till named by Karrow (1963) 
from its type-area in Wentworth County. 

In the southern part of the basin, the Wentworth till appears as a 
buff coloured, sandy to silty sand till (Karrow, 1968), often containing 
stones or boulders. These deposits are largely derived from the parent 
materials of the Lower and Middle Silurian dolomites, shales and sandstones 
that outcrop to the east of the basin. In the northern part of the 
basin, the till becomes reddish in colour due to inclusion of the red 
coloured Ordovician Queenston shale. 

The Wentworth till in this area has largely been deposited as the 


Atal 


Paris, Galt and Moffat moraines, which together form a broad morainic 
belt with a depositional trend from the southwest to the northeast (Map 
2). The Paris moraine lies between the Eramosa River and Blue Springs 
Creek, and the Galt moraine lies to the east of the main Blue Springs 
Creek channel. The Moffat moraine forms the drainage divide on the 
extreme eastern part of the basin. Between these moraines, there are 
glacial spillways which are now occupied by the underfit main channel of 
Blue Springs Creek. The creek near Eden Mills and the main tributary 
south of Acton (Map 1), flow in spillways which have cut through the 
Paris and Galt Moraines, respectively. 


Kame, Esker and Outwash Deposits 


Kame, esker and outwash deposits are glacial meltwater sediments of 
‘stratified sands and gravels. These deposits are often clean washed, 
although they can include till material; they are usually younger than 
the tills and are commonly found as surface deposits. 

Kames are irregular, hummocky accumulations of sands and gravels 
containing some till, and are formed from the sediment carried by 
meltwater on top of the ice. In the Blue Springs Creek basin, kames are 
abundant in the north and west (Map 2). 

Eskers are sinuous ridges of sand and gravel which are formed from 
sediment laden streams flowing under stagnating ice masses. They tend 
to be aligned along the direction of ice movement... Several eskers occur 
in the basin and are shown on Map 2. : 

Outwash plains and terraces are extensive accumulations of sand and 
gravel laid down as relatively level deposits whose surfaces are marked 
by channel scars and kettle holes. In the basin, the outwash sands and 
gravels occur mainly as terrace deposits along the present course of 
Blue Springs Creek (Map 2). 

Kame, esker and outwash deposits have similar lithologies con- 
taining sands and gravels, and are useful sources of building materials. 
An exposed gravel pit in the northern part of the basin shows a strati- 
fied deposit of kame sand and gravel, about 10 feet thick, resting on 
top of the Wentworth till. 

Kame and outwash deposits are usually very permeable and if saturated, 
can provide high yields to water wells. In the Blue Spring Creek basin, 
however, the kame and outwash gravels are shallow surficial deposits and 
are usually above the water table. 


Mechanical and Carbonate Analyses of the Pleistocene Deposits 


Mechanical analyses of samples of the overburden were carried out 
to determine the sand, silt and clay content and to compare these to the 
results obtained by Karrow (1963, 1968) for tills in the Guelph and 
surrounding areas. Nine samples were taken at different depths from 
three observation wells that were being drilled in the area. The 
location of these wells BS-la, BS-2b and BS-3, are shown on Map 3; the 
sample depths, well log descriptions, and the results of the mechanical 
analyses are shown in Table 2. Wells BS-la and BS-3 were drilled 
through surface deposits mapped as kame by Karrow, while the map unit at 
well BS-2b was shown as swamp (Map 2). It was expected that till would 
be recovered at depth, as kame and marsh are commonly surficial features. 

To be consistant with the analyses by Karrow (1968), the mechanical 
analyses figures given in Table 2 are only those for the fraction smaller 
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Figure 4. Mechanical analyses of tills in the Blue Springs Creek 
basin and surrounding areas. 


than 10 mm. Karrow (1974, personal communication) had removed all 
stones larger than about one half inch before performing mechanical 
analyses on the tills collected in the Guelph area. 

The most notable feature of the analyses (Table 2) was the large 
variability of the clay-silt-sand content of samples found in the Blue 
Springs Creek basin. To compare these samples to those found by Karrow 
(1963, 1968) in the Guelph area, and also in the Hamilton and Galt areas 
to the south of the basin, the figures were plotted on a triangular 
plot, Figure 4. Most of the samples from Blue Springs Creek basin have 
less clay than any of the till samples obtained by Karrow, and only 
three samples have compositions similar to those of the tills in the 
Guelph or Hamilton and Galt areas. This shows that most of the samples 
tested in the Blue Springs Creek basin were from thick kame deposits and 
were not till deposits. 

Carbonate analyses were also carried out on the samples from the 
basin and the results are recorded in Table 2. Similar to the mechanical 
analyses, there were larger variations in the carbonate figures from the 
basin samples than in those obtained by Karrow (1968) from the tills in 
the Guelph area. The carbonate analysis figures confirm that the samples 
tested were not till. 


RECENT GEOLOGY 


Recent deposition can be found as alluvial and swamp deposits. The 
alluvial deposits, composed of gravel, sand, silt and clay, are not 
widespread in the basin, occurring only in the Blue Springs Creek valley 
near Eden Mills (Map 2). Bogs and swamps are common in the flat, poorly 
drained areas of the basin and in areas where ground water discharge is 
taking place. Swamps are often covered with shrubs and contain peat, 
muck and marl. 


Soils 


The type of soil formed in an area is related to the type of parent 
material, the drainage and the topography of the area. The soils in the 
Blue Springs Greek basin are derived from the surficial glacial material 
and from the dolomite bedrock exposed near the creek. The soils in the 
area have been mapped by Hoffman et al (1963) and by Gillespie et al 
(1971), and are shown on Map 6. 

The predominant soil is the Dumfries loam, which is derived from 
the Paris and Galt moraine material, and is usually stony and unpro- 
ductive. On the till plains where the topography is more subdued, the 
Guelph loams have developed. These are better soils but can also be 
stoney or may have high relief, hindering crop production. The esker 
and kame deposits and the bedrock exposures near the creek develop a 
variety of poorer soils (Map 6). 

The soils map can be used to delineate moist areas where the water 
table is near ground level. The Dumfries, Killean, Lily, Mesisol and 
marsh soils are a progressively wetter succession of soil types. 

The Mesisol and marsh areas (Map 6) occur in several places in the 
Galt moraine and can be considered as areas where evapotranspiration 1S 
high, even during climatically dry periods in the summer. Marsh and 
Mesisol areas are also considered as possible ground water discharge 
areas. Conversely, well drained soils such as the Dumfries, Or thin 
soils such as the Donnybrook Loam, are indicative of areas of relatively 
low evapotranspiration during dry periods. 


the 


HYDROGEOLOGY 


This chapter presents an analysis of the hydrogeologic factors 
affecting the behaviour of the ground water regime in the basin. The 
data were compiled mainly from water well and observation well records 
maintained by the Ministry of the Environment. Supplementary data, such 
as streamflow and precipitation records, used in arriving at estimates 
of ground water discharges, were obtained from station networks main- 
tained by the Ministry, the University of Guelph and Environment Canada 
(Water Survey of Canada and Atmospheric Environment Service). 


DELINEATION OF HYDROGEOLOGIC UNITS AND PRINCIPAL AQUIFERS 


On the basis of well data, geologic and hydrogeologic cross-sections 
were constructed approximately perpendicular to the axis of the main 
Blue Springs Creek, as shown in Figure 5. Based on these cross-sections, 
three hydrogeologic units were identified, as shown in Table 3. 


TABLE 3. HYDROGEOLOGIC UNITS AND PRINCIPAL AQUIFERS 
IN THE BLUE SPRINGS CREEK BASIN 


Unit Hydrogeologic Unit Princ lpaleaAguLle im 
al kame, Outwash and alluvial 
sand and gravel 

Ssandyetiia: 
reservoir 

2 sandy till 
bedrock 

3 bedrock aquifer 


The kame, outwash and alluvial sand and gravel unit occurs generally 
as lenses or pockets overlying the sandy till unit; however, it forms 
Only minor ground water aquifers in the Blue Springs Creek basin. 

The role of the sandy till unit is that of an aquitard. Sand and 
gravel lenses occur in the sandy till unit, but for the most part are 
not interconnected and have little influence on the hydrogeologic be- 
haviour of this unit. As vertical leakage from unit one to unit two is 
dictated by the permeability of the till, the two units are grouped 
together into one principal aquifer referred to as the sandy till 
reservoir. 

Underlying the sandy till reservoir is the bedrock unit which 
constitutes the major aquifer in the study area. This unit, referred to 
as the bedrock aquifer, is the major source of ground water supply, and 


discharges significant amounts of ground water to the main Blue Springs 
Creek. 


COEFFICIENTS OF TRANSMISSIBILITY AND STORAGE 


In ground water investigations, the coefficient of transmissibility 
(called the transmissivity in this report) and the storage coefficient, 
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are important parameters that control the rate and quantity of the 
ground water flow in a basin. The transmissivity in the Blue Springs 
Creek basin was estimated using various techniques of analysis. These 
were: 

1) Pumping tests on wells (Section 3.2.1) 

2) Specific capacity of wells (Section 3.2.2) 

3) Ground water flow rates (Section 3.2.3) 

4) Ground water modelling (Section 4.4). 

The storage coefficient was estimated by using the following 
methods: 

1) Pumping tests on wells (Section 3.2.1) 

2) Ground water recharge and discharge rates (Section 3.3.1) 

3) Ground water modelling (Section 4.5). 


Hydraulic Properties of the Aquifers Based on Pumping Test Data 


The transmissivities and storage coefficients of the aquifers were 
estimated from two pumping tests carried out in the basin. [In addition, 
the results of two tests carried out near the village of Rockwood, which 
is just outside of the basin, were included. 

There are a number of ways of calculating aquifer constants from 
pumping test data. For this study, three methods of analysis were used, 
all based on the Jacob well equation (Cooper and Jacobs, 1946). The 
methods used were: 

ls) Measurement of drawdown in an observation well during pumping, 

2) Measurement of drawdown of the pumped well during recovery, 

3) Drawdown-distance method, using the drawdowns in observation 

and pumped wells at the end of the pumping period. 


The derivation of the equations and the analysis method are des- 
cribed in Todd (1959) and Wisler and Brater (1959). 

Wells BS-2a and BS-4a were test-pumped in 1966, after being drilled 
for use as observation wells. Wells TW1-67 and TW2-67 were drilled and 
test-pumped in 1967 during a Ministry ground water survey for the 
Village of Rockwood (Sobanski, 1968). The location of these wells are 
shown on Map 3 and the results of the pumping tests are given in Table 
4. 

At TW1-67, there was a marked discontinuity in the drawdown-time 
plot after 20 minutes. This was interpreted as being due to a region of 
high transmissivity in the vicinity of the well surrounded by an area en 
lower transmissivity (Sobanski, 1968). The transmissivities found for 
the bedrock aquifer ranged from 1000 to 12000 IGPD/ft, with a mean of 
approximately 4200 IGPD/ft. As indicated by the range of the transmis- 
sivity figures obtained during the tests, the bedrock aquifer is con- 
sidered to be quite heterogeneous. 

Only one pumping test was carried out in the overburden aquifer. 
The transmissivity found from this test was 6800 IGPD/ft (Table 4). 

The coefficient of storage around the pumped wells was also esti- 
mated (Table 4) and averaged 2x1072. This figure is slightly lower than 
that which is commonly encountered in confined aquifers, where the 
coefficient of storage is usually between 5x107-2 and 500x107° (Todd, 
1959). The coefficient of storage is directly related to the thickness 
and porosity of the aquifer (Murray, 1970); a low storage figure usually 
indicates a low porosity and/or a thin aquifer. As the Amabel Formation 
is fairly thick (approximately 180 ft. near Rockwood) and has a low 
coefficient of storage, it can be considered to have a relatively low 
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porosity. This low porosity is further supported by the low value of 
absorption found for the Amabel dolomite. The absorption is the amount 
of water that a piece of dry rock can absorb, if immersed in water. The 
absorption was measured at 0.93 per cent for the Amabel near Adamsville 
(Hewitt, 1964b). The low porosity and low coefficient of storage 
support the idea of the flow being confined to channels along joints and 
bedding planes, rendering the effective thickness of the aquifer con- 
siderably smaller than the total measured thickness. 


Evaluation of Transmissivity from Specific Capacity Data of Wells 


Transmissivity values were estimated from well records to provide 
an evaluation of the aquifer characteristics throughout the Blue Springs 
Creek basin. The bedrock and overburden aquifers were compared, and the 
transmissivity in the bedrock aquifer was further related to the aquifer 
depths and piezometric surfaces. 

The specific capacity of a well is the rate of discharge of water 
from the well (expressed in gallons per minute), divided by the drawdown 
(in feet). The specific capacity depends on the aquifer character- 
istics, the well diameter, the pumping time, the screening in the well 
and a number of other factors. However, a high specific capacity of a 
well generally indicates a high transmissivity of the aquifer. 

The relation between the specific capacity of a well and the 
transmissivity has been studied by Theis et al (1954). They calculated 
the theoretical specific capacities of wells in relation to the trans- 
missivity and coefficient of storage of the aquifer, the pumping perioc 
and the well diameter, assuming that the well was fully penetrating and 
was fully developed. From the data presented in Theis et al (1963) -the 
relation between the transmissivity and the specific capacity is shown 
in Figure 6, for a well diameter of 0.5 feet, storage coefficients of 
QO. lL and 2x107° and pumping periods of 4 hours and 1 day. 

The transmissivity values near each well were calculated from the 
specific capacity data by using the relationship shown in Figure 6. 

Most of the wells were approximately 6 inches in diameter; if the well 
was considerably larger or smaller, the relationships given in Theis et 
al (1963) were used to estimate the transmissivities. Most of the 
bedrock wells were assumed to be under confined conditions, with the 
storage coefficient taken as 2x10”. Most of the overburden wells were 
assumed to be under water table conditions, with the storage coefficient 
taken as 0.1. Wells having a low specific capacity (less than O.ct 
IGPD/ft) and a short pumping time (less than 2 hours) were excluded from 
the analysis, as the yield from such wells is largely derived from water 
held in storage inside the well itself and not from the aquifer. Due to 
the assumptions inherent in this method, the transmissivity values 
obtained from this analysis should be regarded as approximations. 

In order to determine the median values and the range of values | 
obtained for the transmissivities, a statistical analysis (after Siddiqui 
and Parizek, 1971) was applied to various categories of wells. The 
transmissivity values in each category were listed in ascending order of 
magnitude and assigned probabilities according to the relationship: 
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Figure 6. Relationship between specific capacity and transmissivity. 


100m 
Genin (1) 


where F = percentage of wells where trans- 
missivities are less than the 
Eransmissivity of a well of serial 
number m, 


m = serial number of well arranged in 
ascending order of transmissivities, 


n= total number of wells. 


The transmissivity values were then plotted against percentage o! 
wells on logarithmic probability paper. This technique implies that a 
straight line can be fitted to the plotted points if the observations 
are a random sample from a log-normal distribution. 

The analysis was carried out for wells penetrating the overburden, 
the bedrock, and for bedrock wells which penetrated aquifers at various 
depths and had different static water levels. 

Overburden Wells .... In order to estimate overburden aquifer 
characteristics in the region, data for eight overburden wells in the 
basin and 49 overburden wells from the immediate surrounding areas were 
analysed (see Map 3). 

The logarithmic probability plot for the overburden wells is shown 
in Figure 7, and the number of wells, 10 percentile, median and 90 
percentile transmissivity figures are shown in Table 5. The 10 and 90 
percentile figures are the transmissivities not exceeded by 10% and 90% 
of wells, respectively. These figures give a measure of the dispersion 
of the transmissivity values; a large difference between the 10 and 90 
percentiles indicates a large spread and a high standard deviation. 
standard deviation can be given by (Brookes and Dick, 1951): 


_ *(90) - *(o) 
2, 


oO (2) 


where X,. and X are the 90 and 10 percentile readings, respect- 
ively. thentne ee ee ee are assumed to be log-normally distri- 
buted, the standard deviations which are recorded in Table 5, are of the 
log-transformed data. The log transform of the median figures are also 
recorded in this table for use in statistical analyses (Section 3.2.2.3). 

From Figure 7, it can be seen that the points plot on a straight 
line, supporting the hypothesis that the transmissivity values are 
distributed log-normally. The wide variation in the transmissivity 
figures can be seen from the percentile figures in Table 5, indicating 
the heterogeneity of the overburden aquifer. 

Bedrock wells .... In the basin, 186 bedrock wells were recorded, 
99 of which were selected at random and used for transmissivity an- 
alysis. The transmissivity probability plot for these 99 wells is shown 
in Figure 8 and the characteristics of the plot are recorded in Table 5. 
It can be seen from Figure 8 that the lower transmissivity values plot 
on a straight line, whereas the higher values deviate from this line. 
There are two possible explanations for this deviation. 

The specific capacity of a well is the pumping rate divided by the 
drawdown. At high transmissivities, the drawdown is small and is 
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Figure 7. Transmissivity - probability graph for overburden wells. 
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sometimes reported as zero. In these cases, the specific capacity was 
calculated by assuming the drawdown to be 0.5 feet, introducing a possible 
error. The points where the drawdown was recorded as zero are marked on 
Figure 8; it can be seen that these include most of the high trans- 
missivity figures. 

The other possible explanation is that some of the high trans- 
missivity wells are intrinsically different. The majority of the wells 
are tapping the fissured dolomite; however, a few wells could penetrate 
solutional channels in the dolomite, which would considerably elevate 
the transmissivity values. 

Some of the overburden wells showed zero drawdown during short 
pumping tests (Figure 7) and yet did not display such a deviation at the 
higher transmissivity values. For this reason, the bedrock wells are 
probably from two different populations, with some of the wells tapping 
solutional channels or highly fissured and therefore highly permeable 
portions of the dolomite. 

The eight bedrock transmissivity values obtained from the pumping 
tests (Section 3.2.1) were also plotted on Figure 8. They plot very 
close to the line obtained from the specific capacity of the wells. 

This shows that the transmissivity values obtained from the specific 
capacity data are comparable to the transmissivities obtained .from the 
pumping test data. 

Comparison of overburden and bedrock transmissivity values .... 
From Table 5 it can be seen that the median transmissivity values for 
the bedrock and the overburden wells are similar at 3200 and 3500 IGPD/ft, 
respectively. A Student 't' test (Gregory, 1963) on the log-transformed 
data showed that there was no significant difference between these 
values at the one per cent level. 

The standard deviation of the overburden wells, however, is larger 
than the standard deviation of the bedrock wells (Table 5). The 'F' 
test (Hays, 1973) showed that the difference between the standard 
deviations of the bedrock and overburden wells is significant at the one 
per cent level. 

The median values of transmissivities of the two aquifers are 
similar; the overburden aquifer, however, has more variability than the 
bedrock aquifer. 

Bedrock well analysis .... Further analysis was performed on the 
bedrock wells to determine if the depth where water was found or the 
depth to the static water level had any relation to well yield. The 
bedrock wells were split into four categories according to the depths to static 
water level and depths to water found and plotted in figures 9 and 10. 
The characteristics of the curves are recorded in Table 5. It was found 
that the points in figures 9 and 10) did not: Jie on straight lines, put 
due to the small sample size (generally about 25), it was considered 
likely that the transmissivity values were log-normally distributed, and 
that the deviations seen were random. Standard significance tests were 
carried out on the assumption that the values were log-normally distributed. 

For the depths to static water level, the median transmissivity 
values in each category in Table 5 appear to increase as the water level 
becomes shallower, as was found in Pennsylvanian rocks by Siddiqui and 
Parizek (1971). However, these increases were found to be not Signifi- 
cant at the 1 per cent level. In addition, the standard deviations of 
the figures are not significantly different at the one per cent level 
between any two categories. 
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Figure 8. Transmissivity — probability graph for bedrock wells. 
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The values of depth to water found are plotted in Figure 10 and the 
characteristics of these points are also listed in Table 5. This 
analysis differentiates the wells according to the depths to the water- 
yielding zones in the bedrock. No significant differences were found 
between the median figures; however, a significant difference was found 
between the variability in the categories <69 feet and those >70 feet to 
water found. Thus, the wells tapping shallow aquifers do show a larger 
variation than wells tapping deep aquifers. This finding supports the 
theory of shallow development of solutional channels (Swinnerton, 1932), 
where the range of transmissivity values would be expected to be larger 
than in deeper portions of the rock with less well developed solutional 
channels. It is not known, however, why the median transmissivity 
values in this basin are not significantly different for different 
depths to water found. 


Computation of Transmissivity Values Using Ground Water Flow Rate 


The mean regional transmissivity figure can be estimated by con- 
Sidering the ground water flow rate into the streams. In any aquifer, 
the ground water flow is governed by Darcy's Law, which for discharge to 
streams can be stated as: 


Or=—2TrL (3) 


where Q is the mean gain from ground water flow from 
both sides of the stream along a reach of length L, 


T is the mean regional transmissivity, 


and I is the mean hydraulic gradient towards the 
stream. 


The estimation of the transmissivity was carried out for the section 
of stream between the federal gauge and the Boy Scouts Camp gauge, also 
Shown in Figure 15. This reach is near the centre of the basin, is 
relatively straight and draws its baseflow from the bedrock aquifer. 

The stream length is 2.9 miles and the mean ground water gradient, 
as measured from Map 7, is 84 ft/mile on the west side and 112 ft/mile 
on the east side of the creek. The mean baseflow contribution from 1968 
to 1971 to the creek between the Boy Scouts Camp and the federal gauge, 
was 3.9 MIGPD. 

Using the mean ground water gradient and the measured values for 
the length and baseflow with Equation 3, the mean transmissivity (T) 
was found to be 7000 IGPD/ft. This value is an estimate of the mean 
transmissivity of the bedrock aquifer adjacent to the creek between the 
Boy Scouts Camp and the federal gauge. It can be compared to the median 
values determined in the previous sections (3.2.1 and 3.2.2) on well 
analysis, if the distribution of the transmissivity values is considered. 

The transmissivity values of the bedrock aquifer in this basin were 
found to be approximately log=normally distributed (Section: 3.2.2) sane 
can be shown (Kendall, 1966), that the arithmetic mean of the transmis- 
Sivity is related to the median by the relationship: 
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where xX In’ (T) e5)) 
T is the arithmetic mean of the trans- 
missivity values 


X is the median transmissivity value 


and O. is the standard deviation of the log- 
transformed data 


Using the values for X and 0 derived from Table 5, the mean 
transmissivity of the bedrock aquifer can be estimated, based on the 
specific capacity of the wells as being 7400 IGPD/ft. Thus the arith- 
metic mean transmissivity found from well analysis agrees very closely 
with the transmissivity derived from the ground water flow rate. 


GROUND WATER MOVEMENT 


In most temperate basins containing clastic rocks, the regional 
ground water flow is from the high points in the basin towards the lower 
streams or springs. The streams are usually effluent (Meinzer, 1942) 
and are continuously being recharged by ground water. In any karst 
area, however, the situation can be complicated by the presence of well 
developed solutional channels which can transmit water across topogra- 
phic divides (for example White and Schmidt (1966)) or interfere with 
normal flow of ground water into the streams. 

A piezometric contour map was prepared for the Blue Springs Creek 
area, uSing bedrock water well information, as shown on Map 7. From 
this map the approximate location of the ground water divide was found 
to be not coincident with the topographic divide. The deviation occurs 
in the southeastern part of the basin where ground water is lost from 
the topographic basin, and in the northwestern part of the basin where 
ground water flows into the basin. The ground water basin has an area 
of 26.4 square miles (68.4 km2), which is 13 per cent smaller than the 
topographic basin area (30.2 square miles or 78.2 km) . Because of the 
smaller ground water basin area, the baseflows in the stream will be 
smaller than those calculated using the topographic basin area. 


Monthly Percolation and Ground Water Recharge 


Monthly amounts of percolation and ground water recharge in the 
Blue Springs Creek basin were estimated by determining the soil water 
and the ground water budgets. The hydrologic factors which may affect 
the percolation and ground water recharge include precipitation, evapor- 


ation, snow depth, overland flow, ground water flow, ground water storage 


and soil storage. The interaction of these factors can be illustrated 
in the simplified conceptual model shown in Figure 11. In this figure, 
there are four storage zones (the snowpack, the soil, the lower soil and 
the ground water), where the conservation of mass equation can be 


applied: 


ao 


PRECIPITATION (P) 


TRANSPIRATION (E) 


SNOWPACK 
STORAGE (SP) 


SNOWMELT 


WATER AVAILABLE 
ABOVE SOIL (WwW) 


SOIL OVERLAND 
STORAGE (SS) FLOW (0) 


PERCOLATION (PC) 


LOWER SOIL 
STORAGE (LSS) 


RECHARGE (R) 


GROUND WATER BASEFLOWS (B) 
STORAGE (GWS) 


STREAMFLOW (S) 


Figure Il. Conceptual model used to estimate percolation and ground water 
recharge. 
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in the snowpack 


P = W + ASP; (6) 
in the soil 


W= PC +0+E + ASS; (7) 
in the lower soil, | 
FC_£E= R + ALSS; (8) 
and in the ground water 
R = B + AGws (9) 


where P is the precipitation, 
W is the water available above the soil, 
ASP is the change in snowpack storage, 
PC is the percolation below the soil zone, 
O is the overland flow, 
E is the evapotranspiration, 
ASS is the change in soil storage, 
ALSS is the change in lower soil storage, 
R is the ground water recharge, 
B is the baseflow, 

and AGWS is the change in ground water storage. 


The evapotranspiration in Equation 7 implicity includes any snow- 
pack and ground water evaporation. 

The mean monthly recharge was estimated for the years 1966 to 1972, 
by solving equations 6 to 9. Much of the data in these equations are 
available from existing records; the precipitation was taken from the 
meteorological records, the evapotranspiration was calculated by the 
Thornthwaite (1948) method, the snowpack storage was taken as one tenth 
of the snowpack depth at the Blue Springs Creek basin, and the stream- 
flow at the federal station was separated into the baseflow and overland 
flow components. These data are listed in Appendix 1. Using these 
data, the only parameters in equations 6 to 9 that are unknown are the 
soil storage, the percolation, the ground water storage, the lower soil 
storage and the recharge. 

The soil storage was estimated for each month by making the follow- 
ing assumptions concerning the percolation and maximum soil storage 
levels. The soil was allocated a maximum storage of the field capacity, 
and any excess water available was allowed to percolate. If the soil 
was below field capacity, the percolation was considered to be zero. 

Twenty measurements were made of the moisture storage in the top 
four feet of the soil zone at three sites during 1969 in the Blue 
Springs Creek basin. The measured values range from about 10 inches in 
the spring to three inches in the late summer. It was found that these 
soil moisture measurements were closely modelled if the field capacity 
was set at 10 inches. Using this figure for the field capacity, the 
monthly percolations were determined from Equation 7. The data are 
shown in Appendix 1 and the mean annual values are shown in Table 6. 

The monthly percolations are plotted in Figure 12. 
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The percolation estimates obtained by this analysis must be con- 
sidered to be approximate values, due to the number of assumptions made 
using this method. In particular, the model did not represent conditions 
in the summer very well. The model predicted that the percolation in 
the summer was generally zero, and during the summer the mean monthly 
evapotranspiration was usually larger than the precipitation. However, 
a heavy rainfall during the summer would probably give some percolation. 
To model these short events would require reducing the time step to 
about a day. 

The recharge can be estimated by using the ground water balance 
(Equation 9) , if the ground water storage can be determined. The 
change in ground water storage was approximated as: 


Acws = AWL.SC.A (10) 


where AGWS is the change in ground 
water storage, 


A WL is the change in mean 
water level in the area, 


SC is the mean coefficient of storage, 
and A is the basin area. 
Combining equation 10 and 9: 

R= B + A(WL).SC.A Ci) 


The mean water level in the basin was estimated by taking the mean 
level in the six observation wells representing discharge and recharge 
areas in the basin. 

During some of the summer the recharge to the ground water is zero. 
During these periods the hydrographs of the observation wells can be 
seen to decline approximately exponentially, without any significant 
rise in water level caused by recharge. By inspection of the observation 
well records, 14 summer months (Figure 13) were identified when the 
recharge appeared to be zero. 


If the recharge is zero, Equation (11) can be restated as: 


SVE! ee -B 
AWL.A 


(¥2} 


Using Equation 12 for the 14 summer months when the recharge was 
apparently zero, the value of the coefficient of storage was found to 
vary between 0.03 to 0.05, with a mean of about 0.04. This mean value 
of storage coefficient was then used in Equation 11 to get the monthly 
recharge values for the years 1966 to 1972. These values have been 
plotted on Figure 13 and are listed in Appendix 1. For a few months the 
recharge was calculated to be negative; these values were set to zero, 
and the negative value was subtracted from the following month's recharge. 
For several of the summer months when the recharge was considered to be 
zero from inspection of the well hydrographs, a small recharge value was 
obtained by using Equation 1l. These were due to the approximations 
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inherent in the method, and for consistency these figures were left at 
the value calculated by using Equation ll. 

Comparing the percolation (Figure 12) with the recharge (Figure 
13), it can be seen that the percolation is considerably more variable 
than the recharge. In addition, the percolation peaks are about one 
month ahead of the high recharge figures. The lower soil zone, there- 
fore, acts to smooth and delay the percolation before it reaches the 
ground water. 


Ground Water Discharge Zones 


Ground water discharge zones were investigated by three different 
methods. Discharge and recharge zones along the stream were found by 
streamflow measurement. Other discharge zones were found by examining 
the data from the three piezometers and the general piezometric map. 

The methods and results are discussed in the following subsections. 

Discrete recharge and discharge zones .... In temperate basins 
streams are usually gaining, but in karst basins the streams often 
recharge aquifers and receive ground water discharge at different points 
along their course. To estimate the amount of recharge and discharge in 
different reaches of Blue Springs Creek, the streamflow was measured at 
a number of points along the channel, on December 13, 1973; at this time 
flows were relatively low and receding. These measurements are shown in 
Figure 14 (in pocket), both in cfs and in inches of runoff per year. To 
compare the reaches of the stream that may be losing water to aquifers 
and those that gain ground water, a theoretical streamflow was calculated 
for each measuring point by assuming that the runoff was constant over 
the whole basin and was equivalent to 10.5 inches per year as measured 
at the downstream station. If runoff were evenly distributed all reaches 
would be supplying the theoretical streamflow. The differences between 
measured and theoretical values were used as an indicator of recharge 
and discharge areas and are shown in Figure 14. 

The major area of losing flow appears to be the reach in the 4th 
Concession, Nassagaweya Twp., where there is a difference of approxi- 
mately 4.3 cfs between the actual and theoretical streamflow. This 
water enters the ground water and it can be seen from Figure 14 that a 
comparable quantity of water discharges into Blue Springs Creek farther 
downstream. Other reaches of low runoff can be seen in the southerly 
part of the basin. These may be due to local infiltration to underlying 
aquifers or to ground water flowing to more distant reaches of the 
streams. All other reaches were found to be gaining, with runoff 
greater than the theoretical runoff. 

Field surveys were carried out to check for sinks and resurgences 
along Blue Springs Creek and throughout the basin. A number of these 
sites were identified and are shown in Figure 20. Several streams sink 
in sand and gravel along the streambed, although no discrete sink points 
were seen that carried water all year. 

There are a number of springs in the basin (Figure 20), most of 
which are small. One significant spring in the area can be found on the 
Boy Scouts Camp property, as the Blue Spring. This spring was discharg- 
ing 0.47 cfs on December 13, 1973, which would correspond to a runoff 
from approximately 0.6 square miles, using the runoff rate at Eden 
Mills. The topographic basin area contributing to this spring is only 
about 0.05 square miles, and from the piezometric map (Map 7), it can be 


seen that the source of this spring probably includes the higher land to 
the southeast. 
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Recharge and discharge areas from well piezometer data .... The 
recharge and discharge areas were also investigated near the piezometer 
nests in the basin. A piezometer nest consists of a well containing a 
number of tubes which measure the hydraulic head in aquifers at dif- 
ferent depths. From the static levels in each tube, the vertical 
component of ground water movement can be inferred. Three piezometer 
nests are located in the basin and are shown in Figure 15 (in pocket) as 
wells BS-lb, 4b, and 5 containing 4, 2 and 5 piezometer tubes respect- 
ively. The static water level measurements for the wells during 1972 
aresplotted in figures 16, 17 and 18. From these graphs, the average 
level in each piezometer tube in 1972 was calculated and converted to 
relative static levels, which are shown in Figure 19. The inferred 
vertical flow directions, based on the mean static level in the pie- 
zometers, are also shown in Figure 19. 

In well 1b, located in the northern part of the basin, the pre- 
dominant flow direction is upwards, indicating a discharge zone. The 
discharge of ground water causes a marsh to be formed in the vicinity of 
the well. This well, however, is near both the topographic and ground 
water divides (Figure 15) where, under normal conditions, the ground 
water would be expected to be flowing downwards. The lower aquifers 
near this well must be recharged from areas of higher elevation, probably 
from the high region to the north of the basin. Further information on 
the type of recharge can be obtained by inspection of the piezometer 
levels, Figure 16. During dry periods when the water levels are falling, 
the four piezometer tubes indicate that the water is flowing upwards and 
discharging at the surface. During wet periods, however, when the water 
levels are rising, the second deepest piezometer tube tapping a gravel 
bed at the 128-140 foot level, has the highest head. This bed must be 
connected to the recharge area by a relatively permeable path, while the 
aquifers at other levels must be more localized, and have a poorer 
connection to the recharge areas. 

In well number 5, the upper piezometer tube had been filled with 
gravel to the 9.4 foot depth, and during much of the year no readings 
could be taken as the water level was below 9.4 feet. 

However, when the static level is high enough for readings to be 
taken, the level is always at least three feet higher than the level in 
the other piezometers and this difference was assumed to exist through- 
out the year. From Figure 19, the flow direction is predominantly 
downwards, indicating a recharge zone. 

Well number 4b, is also in a recharge area. The flow direction, as 
shown in Figure 19, is downwards. 

Delineation of discharge zones on the basis of water level informa” 
tion .... In order to determine the approximate boundaries between the 
regional discharge and recharge zones from the bedrock aquifer, the 
following information was used. 

A) The piezometric surface map (Map 7) was compared to the topo- 
graphic elevations. The areas where the piezometric surface 
was at or above the ground surface indicate discharge areas; 

B) The topographic map (Map 1) was inspected for marshy areas; . 

C) The soils map (Map 6) was inspected for marsh and Mesisol soil 
types. These areas occur either where ground water is being 
discharged or in depressions where the drainage is poor; 
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Figure 20. Approximate discharge and recharge zones in the Blue 
Springs Creek basin. 
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D) The air photographs of the region were studied for moist 
areas. These marshy areas were found to be darker on the 
aerial photographs and could be distinguished from the lighter, 
cultivated areas; 

E) The piezometric well data (Section 3.3.2.2) were used; 

F) In the field, springs were identified and located. Much of 
the field locating was carried out by Dickinson and Whiteley 
(19°70) 

G) The chemistry of the ground water was used to estimate re- 
charge and discharge zones. This method is described in 
Sections... 

On the basis of the above evidence, the discharge and recharge 
zones were estimated and have been marked on Figure 20. The area of the 
discharge zones was measured at 1.9 sq. miles, (4.9 sq. km) or about six 
per cent of the basin area. 


STREAMFLOW AND BASEFLOW CONDITIONS 


The mean annual streamflows for the calendar years 1967 to 1971 are 
given in Table 7 for three stations in the Blue Springs Creek basin. 

The stations are the Cedar Grove and Boy Scouts Camp stations maintained 
by the University of Guelph, and the federal gauge 2GA-31. The station 
locations are shown in Figure 15. The topographic drainage areas above 
each station, the runoff per year, the annual precipitation and the 
percentage of runoff to precipitation are listed in Table 7. 

The runoff from the federal gauge and the Boy Scouts Camp gauge is 
approximately 40 per cent of the precipitation. The runoff from the 
Cedar Grove gauge however, is only 23 per cent of the precipitation; 
this lower figure is due to the relatively smaller amount of ground 
water discharge in the Cedar Grove basin. 


Ground Water Discharge Based on Hydrograph Separation 
In any hydrograph, such as the one illustrated in Figure 21, it is 


generally considered that there are four major components of the stream- 
flow (Linsley et al, 1949). These are: 


a) Channel precipitation, which falls directly on the stream 
channel. 
b) Surface runoff, which is conventionally regarded as the water 


that flows over the surface of the land. ‘This-is the major 
component of flood flows. 

fei Interflow, which is the water that flows laterally through the 

soil without reaching the water table. 

d) Baseflow, which is derived from the ground water reservoir. 

This is the major component of dry weather flows. 

In practice, the separation of the hydrograph is difficult as there 
are no strict divisions between the components. Thus, some water may 
flow on the surface initially, then in the soil for a short distance, 
and back over the surface before reaching the stream network. This 
water may be classified either as interflow or surface runoff. 

For this study, the streamflow was separated into the baseflow and 
the sum of the other three components, which together were called the 
overland flow. This analysis was carried out by conventional baseflow 
separation techniques (Gray and Wigham, 1970). 

The hydrograph of streamflow, during periods when all discharge is 
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derived from ground water sources, is known as the baseflow recession 
curve. Several theoretical equations have been proposed to describe the 
recession curve. They are generally of the form (Barnes, 1940) : 


Q. = OR K =-t ES) 
where: 
Q, = discharge at any instant, 
Qo = discharge at some initial time, 
K = recession constant, 
and t = time interval between On and Or 


A baseflow recession curve within the period extending from 
February 15 to March 12, in the year 1969, at the federal gauging 
station, was selected and assumed to represent the ground water outflow 
in the basin, as shown in Figure 24. A winter recession curve was 
selected because evapotranspiration losses to the atmosphere during that 
period were assumed negligible. Substituting the flow values from the 
recession curve from Figure 24 into Equation 13 gives a recession con- 
stant K of 0.989, or a 'half life' for the recession of 64 days. Thus 
the baseflow after a 64-day period will be reduced by one half if the 
aquifer has not been recharged. 

The slope of the recession curve was drawn on the stream hydrographs 
(Figures 22 to 24) showing total runoff on semi-logarithmic paper. The 
recession curve plotted in this manner appeared to fit the falling 
stages of the hydrographs for the time periods 1966-1969. In order to 
determine the limits of the baseflow recession curves on the stream 
hydrographs, observation well (2b) was used as an index well which was 
assumed to reflect the natural dewatering and recharge effects on the 
ground water reservoir. 

Figures 22 to 24 illustrate the manner in which the stream hydro- 
graph was separated into two components, surface runoff and ground water 
discharge. The daily ground water discharges beneath the recession 
curves were then added to give the total monthly ground water discharge; 
Table 8 summarizes the monthly ground water discharge values. 

The lowest and highest monthly ground water discharge, based on the 
three year period, were 8.5 and 25.6 cfs, respectively. The mean ground 
water discharge was 14.0 cfs, which was 68% of the total streamflow. 


Specific Yields from Different Sub-basins in the Blue Springs Creek Basin 


When an aquifer is drained under the action of gravity, the amount 
of water yielded is less than the total volume of the aquifer material 
drained (Wisler and Brater, 1959). The quantity of water that is yielded 
from the aquifer is called the specific yield, and is usually expressed 
as a percentage of the total aquifer volume. During dry periods when © 
the aquifer is not being recharged, and if evapotranspiration and pump- 
age from the aquifer are small, the specific yield can be given as: 
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TABLE 8. MONTHLY AND MEAN MONTHLY GROUND WATER 
DISCHARGE (IN CFS-DAYS) AT THE FEDERAL GAUGE 
02GA031 IN THE BLUE SPRINGS CREEK BASIN 


ee! ee eS ee ee eee ee 
Water Years 


Month LI6G=1IG LOO7>06 LIC6S=—69 Mean 


October LSOGs BWAa eg! Sidisn 7 20004 
November ei Day 4) Lg S Ons Do 2a 
December Sa Se ADU I 5 Os0 437.6 
January SHOLO) Fi A205 569.3 430.1 
February 366.6 5S LSs5 596.7 493.9 
March Soonc 695.2 586.5 BSS) 5 Z 
April 693.2 aT Ed ie SG dal ote. 
May 5OSe2 507-5 Pe eseee 5S /Ge2 
June AVGS9 S336 5 480.1 AL iLILs (3 
July 466.2 242.9 349.6 B52 6S 
August Sy Bite! 23555 ZIOn 216 
September 2 ee SoO= 7 18554 INGO te: 
TOTAL Al 3\'5)'5) 5 1 Les. 0 5849-9 SPAS T 8} 
MEAN 

MONTHLY 362.9 AZ2EO0 487.4 BOT 3 & 
Mean 

Baseflow 

Gis Leg Age 1S) AO) 
Mean Stream- 

flow at 

O2GA31 cfs Lows 2Oo8 DA KS DORIS 


Percentage of 
Baseflow to 
Streamflow 64% 68% WALES 68% 


Dee eiis 


a cae x 100 (14) 


where y is the specific yield in per cent, 


Q is the mean baseflow discharge during 
the period At, 


Ah is the mean change in piezometric 
head over the area A. 


The specific yield can be calculated for aquifers in a sub-basin if 
the mean decline in head and mean baseflow can be estimated. 

The specific yield was estimated for the six sub-basins marked on 
Figure 15. These are the Cedar Grove, Boy Scouts Camp, federal gauge 
and the Blue Springs Creek basins above each gauge, and the Knatchbull 
and the Eden Mills sub-basins between the lower two gauges. 

The calculations were carried out for the period February 15 to 
March 12, 1969 when the precipitation and evapotranspiration were assumed 
to be negligible so that the entire streamflow was taken to be baseflow. 
The change in head in the aquifer was estimated from the recording 
observation wells la, 2a, 2b, 3 and 4a and at the piezometer nest 5 (see 
locations on Figure 15). The heads in these six wells are shown in 
Figure 25. The mean change in head in any sub-basin was calculated by 
taking the average decline in head of the wells in or near that sub- 
basin. The wells considered for each sub-basin and the average decline 
in head are shown in Table 9. 

The streamflows at the four gauging stations are plotted on Figure 
26. The mean baseflows were calculated by averaging the daily flows. 
The baseflows in the Knatchbull and Eden Mills sub-basins were taken as 
the differences between the baseflow at the gauge downstream and the 
gauge upstream of the sub-basin. The specific yields were calculated by 
using Equation 14 for the six sub-basins, and are recorded in Table 9. 

The specific yields in the sub-basins range from 7.2 to 12.7 per 
cent (Table 9). The lowest specific yields are from the Boy Scouts Camp 
and Cedar Grove sub-basins, which are in the northern portion of the 
area. The major portion of the ground water in these sub-basins is 
derived from the sandy till reservoir. The highest specific yields are 
found in the Knatchbull and Eden Mills sub-basins, which derive much of 
their ground water from the dolomite aquifer. However, as much of the 
dolomite aquifer is overlain by saturated till, the specific yields 
found in these sub-basins also represent the conditions in the overlying 
till. The sandy till aquifer is composed of kame, esker, outwash and 
alluvial sands and gravels and sandy till. The sands and gravels 
generally have a higher permeability than tills and could be expected 
touhavesa nigher specific yield than the finer till (Todd, 1959), if 
conditions of thickness and degree of saturation were similar. 

The proportion of the surficial sands and gravels (shown on Map 2) 
in each sub-basin is given in Table 9. With the exception of the Cedar 
Grove and Eden Mills sub-basins, the specific yield generally increases 
as the proportion of sands and gravels in the basin increases. In the 
Cedar Grove sub-basin, which has a small specific yield, the proportion 
of surficial sands and gravels is high, but these sands and gravels tend 
to be thin and are not very persistent. The shallowness of these deposits 
can be seen in the geological cross sections A-A' and B-B' of Fig. 5 (in 
pocket). On the other hand the Eden Mills sub-basin, which has a high 
specific yield, contains a small proportion of sands and gravels but 
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Figure 25a. Storage conditions in the sandy till reservoir as shown by 
Observation wells la, 2q,2b in the Blue Springs Creek 


basin, February 15 to March l2, 1969. 
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Figure 25b. Storage conditions in the sandy till reservoir as shown by 
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observation welis 3, 4a and at the piezometer nest 5 in the 
Blue Springs Creek basin, February I5 to March 12, 1969. 
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Figure 26. Daily discharges at the Cedar Grove, Boy Scouts Camp, 


federal and Eden Mills gauges for the period February |5 
to March) 12, 1969). 
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deposits are thick and persistent, as can be seen by cross-section G-G'. 
Therefore, it appears that large specific yields in a basin are related 
to large, thick and persistent areas of kame, esker, outwash and alluvial 
sands and gravels. 
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GROUND WATER MODELLING 


As described in previous chapters, the transmissivities, coeffi- 
cients of storage and recharge of an aquifer are important in under- 
standing its flow pattern and defining its potential yield. In this 
chapter, two phases of ground water modelling are applied to evaluate 
and extend the results obtained by conventional techniques, and to 
provide a better understanding of the regional distribution of these 
aquifer characteristics and their interaction with the streams in the 
basin. 
In the previous chapter, estimations of the transmissivity values 
were obtained, based on well record analysis and on ground water flow 
rates. Although these techniques yielded fairly consistent estimates of 
the transmissivity, they could not be used to obtain satisfactory figures 
for the regional transmissivity throughout the basin. Well record 
analysis provided values for the local conditions near the wells. Analysis 
of the ground water flow rate did provide regional transmissivity values; 
however, these were limited to only parts of the basin, due to lack of 
extensive data. One method of further deriving the regional variations 
of transmissivity is through the use of modelling. In the Blue Springs 
Creek basin, a steady state, finite difference model was used to estimate 
the regional transmissivity values; the figures obtained were then 
compared to those found in Chapter 3. 

The coefficients of storage and recharge to the ground water 
aquifer are other important parameters in understanding the ground water 
regime. Both parameters, however, are difficult to determine by conven- 
tional techniques. A dynamic finite difference model was used to estimate 
the coefficient of storage throughout the basin and to verify the time- 
varying recharge rates to the aquifer calculated by water balance 
considerations. 

One major concern in applying ground water models to simulate 
limestone or dolomite aquifers is the possible presence of karst develop- 
ment and the need to incorporate the associated flow characteristics 
into a modelling routine. As illustrated previously, if the karst 
conditions in a basin are well developed, solutional channels may transmit 
large quantities of water through the aquifer. If such an area is to be 
modelled accurately, the individual solutional channels have to be 
identified and located in the field in order for the simulation model to 
adequately account for their effect. If these channels, however, were 
not identified, the ground water model would not adequately simulate the 
aquifer characteristics. In the case of the Blue Springs Creek basin, 
the karst development has been shown to be of limited extent, the 
effects of the minor solution channelling, sinkholes and springs having 
no detectable effect on the piezometric water level surface of the area 
(Section 2.1.4). For the Blue Springs Creek basin, the ground water was 
modelled by assuming that the flows obeyed the Darcian flow equations. 
This appears valid, considering that the transmission of flow through 
solutional channels is minimal and only local. The model then was able 
to simulate the average aquifer characteristics on a distributed area 
basis, showing their general distribution and variation throughout the 
basin and its surrounding area. 


DESCRIPTION OF THE MODELS 


The programs used to solve the ground water flow equations by the 


aM) 


finite difference method are described by Smith, 1965, and Prickett and 
Lonnquist, 1971. A brief summary of the mathematics and the programs is 


presented here. . 
The generalized two dimensional differential equation of flow in an 


aquifer is given by: 


d dh | dh dh 

—- —) +— — +I=s— LS) 

dy es ay dx ee: ae ac ( 
where T , T are the transmissivity components in the 


x*and y directions respectively, 


I is the sink or source term, representing 
pumpage or infiltration, 


h is the piezometric head in the aquifer, 
S is the storage coefficient, 
and Ceicetime. 


In the steady state case, when the head does not vary with time, 
the right hand side of the equation will be zero. 

Equation 15 is the basic equation for all two dimensional ground 
water flow problems. There is no general solution for this equation, 
but approximate solutions can be obtained by a number of methods. In 
the finite difference method of solution, the differential terms in 
Equation 15 are replaced by their finite difference approximations. The 
aquifer being modelled is then discretized into a nodal network and the 
finite difference equations are solved at each node in turn. 

In order to solve the problem more expediently, some physical 
idealizations have to be made (Figure 27). The aquifer is assumed to be 
thin and confined below by an aquitard, so that the flow in the aquifer 
is assumed to be horizontal. Infiltration can occur to the aquifer 
through the leaky confining layer above it and pumpage or spring dis- 
charge can occur at any node in the model. Recharge or discharge can 
occur through the permeable beds of the rivers which are assumed to have 
a constant head at all times. In practice, these assumptions appear not 
unduly restrictive. 

The first step in the modelling process is to superimpose a grid 
network over the area to be modelled. The grid used for the Blue Springs 
Creek basin and the surrounding area is shown in Figure 28; it was set 
up unequally spaced with more nodes in the basin area where detailed 
head information is available. 

At each node, all sources of water that can enter or leave the node 
are considered, as shown diagrammatically in Figure 29. There are seven 
routes that the water can enter or leave any central node (marked 'O' on 
Figure 29). These are: to or from the four adjacent nodes, to or from 
the river, to or from storage and by direct pumpage or recharge. By 
altering the theoretical head at the central node, it is possible to 
balance the amount of water entering this central node with the amount 
leaving. By repeating this process for every node in the basin, on an 
iterative basis, the theoretical head distribution throughout the aquifer 
can be calculated. 


In the steady state case, the model uses the "overrelaxation" 
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Figure 27. The idealized aquifer system modelled. 
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Figure 28. Node network and transmissivity values used for the ground water model. 
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6.1 


method of solution for the equations (Southwell, 1946) where the head at 
the central node ("O" in Figure 29) is overcorrected according to the 
heads at the adjacent nodes. For the dynamic case, the solution is more 
efficiently found by using an alternating direction implicit method of 
solution (Peaceman and Rockford, 1955), where the heads along any row or 
column are found by consideration of the heads at adjacent rows or 
columns along with the pumping and storage stresses. These techniques 
have been written into two computer programs that serve to solve the 
steady state and dynamic ground water problems. (Pricket and Lonnquist, 
BOTA). 


Model Testing 


A number of tests were carried out in the models to verify the 
programs. A number of runs were made with both the steady state and the 
dynamic model on simple well problems, where the drawdowns could be 
calculated theoretically. The programs were found to model these condi- 
tions accurately. 

Using the initial values for the Blue Springs Creek data, a sensi- 
tivity analysis was also carried out. The effect on the output of a 
change in any of the input values was tested in a qualitive manner. 
Table 10 shows the effect of an increase in the transmissivity, recharge, 
storage coefficient or river permeability on the output of the model. 
The results from these tests were used during the later optimization 
runs for the Blue Springs Creek basin. 


PREPARATION FOR THE MODEL 


To perform steady state modelling, the following data are required 
at each node: 


a) transmissivity (IGPD/ft), 
b) recharge (IGPD/£t2) or pumpage (IGPD), 
CG) river head elevation (ft) and the permeability between 
the river bottom and the aquifer (IGPD/f£t*), 
d) heads (ft) at each node for verification of the model. 
For the dynamic model, further data are required: 
e) storage coefficient (Imp. gal/ft3), and the recharge rate each 


month, and the heads at observation wells for each month. 

The model was run using the above data and the simulated and the 
measured heads were compared. On the basis of any discrepancies between 
the known and modelled heads, the input data to the model were adjusted. 
Initial estimates for the input data were obtained from the hydrogeo- 
logical calculations carried out in Chapter 3; their uses in the model 
development are described below. 
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TABLE 10. RELATIVE EFFECTS OF A CHANGE IN THE INPUT 
PARAMETERS ON THE MODEL OUTPUT 


Increase in 
values of: 


Resultant Change in Output: 


Steady State 
Baseflow 


Baseflow 
Fluctuations 


Steady State Water Level 
Water Level Fluctuations 


Transmissivity None None Reduced Reduced 

Recharge Increased Increased Increased Increased 
Storage 

Coefficient None Reduced None Reduced 


River 
Permeabilit 


None Reduced Reduced Reduced 


Transmissivity 


The transmissivity data were derived from the specific capacity 
data of the wells in the area (Section 3.2). The aquifer was found to 
be very heterogeneous and wide variations between the transmissivities 
of nearby wells were often found. 

When considering regional ground water flow patterns, the mean 
transmissivities of different areas in the basin are important. Thus 
the modelled area around the Blue Springs Creek basin was divided into 
several large rectangular areas (Figure 28) and the mean transmissivity 
in each area was calculated from the wells in that area for initial 
input to the model. 

The presence of solutional channels in a karst area can upset the 
calculated transmissivities and the flow patterns. In the Blue Springs 
Creek basin the karst was not considered to be well enough developed to 
seriously perturb the transmissivity values used in the model. 


Recharge to the Ground Water 


Estimates of the ground water recharge were determined in two ways: 
from considerations of the mean baseflow to determine the long term mean 
recharge for input to the steady state model, and by examining the water 
balance in the basin to determine the time varying recharge rates for 
input to the dynamic modelling. 

The mean baseflow from the Blue Springs Creek basin above the 
federal gauge, for the three years 1966 to 1969, was determined to be 
14.0 cfs (Section 3.4.1). The mean recharge was assumed to be equal to 
the long term mean baseflow; this figure was used as the steady state 
model input. 

The effect on the aquifer of pumpage from wells was determined. 
There are about 200 recorded water wells in the Blue Springs Creek 
basin, all of which are farm or domestic wells, many of which are now 
abandoned. Assuming that all of these wells were pumped at an average 
of 200 IGPD, the total pumpage would be less than four per cent of the 
baseflow determined above. For this reason the pumpage effect of the 
wells on the aquifer system was considered to be negligible. 

The time varying recharge rates were calculated by considering the 
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water budget in the basin, as described in Section 3.3% Leese Ln eis 
section the recharge was calculated by considering the ground water 
budget. The monthly recharge figures used are shown in Figure 13. 


River Heads and Permeabilities 


In the Blue Springs Creek basin all of the water entering the 
ground water was assumed to ultimately discharge into the streams as 
baseflow. For the model input, stream elevations and the permeability 
between the streams and the aquifer were determined. 

The elevations of the streams in the basin and the surrounding area 
were determined from topographic maps. As the stream course rarely 
coincided with any node in the superimposed model grid network, the 
stream system had to be "squared up" to run through the nodes, as shown 
in Figure 30. The stream elevations at the points nearest to the river 
nodes were then determined from the maps. 

The streambed leakage factor between the river and the aquifer is 
difficult to determine. The streambed leakage factor at a node can be 
given as: 


_ Plw (16) 
SL = a 
where SL is the streambed leakage factor at the node, 


P is the mean permeability of the river bed 
material, 

1 is the length of river associated with the 
node, 

w is the mean river width, 

m is the average thickness of the river bed 
material, 

and A is the area associated with the node. 


The thickness and permeability of the river bed material are dif- 
ficult to determine. Therefore, Equation 16 could only be used to 
derive a general index for the streambed leakage factor. After some 
tests with the model, an effective permeability of 2 IGPD/f£t* was assigned 
rOnthe rivers, 1.0 IGPD/£t2 to large creeks, 0.5 IGPD/ft2 to medium 
creeks and 0.1 IGPD/ft* to small creeks. These permeabilities represent 
very good hydraulic conductivity between the streams and the aquifer, 
and effectively render the heads in the aquifer at the stream nodes 
similar to the elevations in the streams. 


Piezometric Heads 


The heads in the aquifer, which were required for model verification, 
were obtained from the water well information. The long term mean heads 
at each node inside the basin were derived from the piezometric surface 
map (Map 7). Data on the change in elevation of the water at six obser- 


vation wells were extracted on a monthly basis and used for verification 
of the dynamic model. 


Storage Coefficients 


The storage coefficient was estimated at about Pea from pumping 
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Figure 30. Modelled stream system heads and permeabilities 
in the Blue Springs Creek basin. 
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tests on wells near Rockwood (Section 3.2.1). This figure tends to 
indicate confined conditions. The storage coefficient was also esti- 
mated at 4x1072 from water balance considerations explained in Section 
3.4. This latter estimation indicates water table conditions. The 
model was run using both estimates for the storage coefficient. 


STEADY STATE MODEL RESULTS 


During initial runs of the model, it was found that the variations 
in the heads over the basin were larger than the measured observed 
variations. This difference could be corrected by increasing the trans- 
missivity values used in the model, or by decreasing the recharge rate 
over the whole basin (Table 10). As the baseflow from the model is 
proportional to the infiltration rate, and was found to be approximately 
correct during the initial runs, the variations in head were conse- 
quently subdued by increasing the model transmissivities. It was found 
that increasing the transmissivities by 1.5 times their calculated 
values subdued the heads sufficiently over a large part of the basin. 
The higher transmissivities in the model, compared to the values found 
from the wells, are concluded to be due to the nature of the aquifer. 

In many consolidated rocks the majority of the flow takes place in 
joints, fissures and bedding planes. Many wells may penetrate poorly 
fissured portions of the rock and so record a lower transmissivity value 
than the true regional value. It appears that the regional transmis- 
sivity values in Blue Springs Creek basin are therefore about 1.5 times 
the median values found from the wells in the area. 

After adjusting the transmissivity values in this way, errors were 
still found in the modelled heads near Blue Springs Creek, where the 
heads were too subdued, and in the northern and southern part of the 
basin, where the modelled heads were too high. These discrepancies 
could be corrected by adjusting either or both of the modelled transmis- 
Sivities or recharge rates near the creek and in the southern and 
northern portions of the basin. It was found that the baseflow at the 
Cedar Grove gauge (located in Figure 15), which drains the northern part 
of the basin, was smaller (5.3 inches a year) than the baseflow at the 
federal gauge (9.8 inches a year), which drains the central and northern 
part of the basin, indicating that the discharge rate is less in the 
Cedar Grove sub-basin than in the basin as a whole. The errors in the 
modelled heads were therefore corrected by adjusting the recharge rates, 
while leaving the transmissivities at their previous values. As the 
Cedar Grove sub-basin generally has at least 50 feet of overburden above 
the dolomite, the recharge was decreased to 7 inches per year in the 
regions having more than 50 feet of overburden above the aquifer (Map 
5). To maintain the measured baseflow from the federal gauge sub-basin, 
the recharge was increased to 15 inches per year over the exposed dolomite, 
where the infiltration capacity is higher than in regions covered by 
thick till. The final recharge rates used are shown in Figure 31. 

Although adjusting the recharge rates over the basin gave a satis- 
factory head distribution from the model, consideration of the rela- 
tively low runoff rate from the Cedar Grove sub-basin (Table 7) suggests 
that there may be significant ground water flow from the northern to the 
southerly parts of the basin. Such a flow is consistent with the 
piezometric information, and would be possible in karst areas, but 
further work would be needed to prove that such a flow exists. This 
could be modelled by including a large anisotropic transmissivity in the 
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Figure 3I. Long term mean effective recharge rates and location of the observation 
well nodes. 
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north-south direction, to allow significant flow to occur from north to 
south. However, such a model was not tried, as adjustment of the recharge 
rate over the basin gave satisfactory results. 

Using the adjusted transmissivity and recharge figures, the modelled 
heads and a histogram of the differences between the modelled and measured 
heads for the nodes inside the basin, are shown in Figure 32. At the 
majority of the nodes (approximately 70 per cent), the differences were 
less than 10 feet, and the largest difference was 28 feet. The modelled 
average baseflow in the federal gauge sub-basin was 9.79 inches per 
year, compared to the measured baseflow of 9.8 inches per year. 


DYNAMIC MODEL RESULTS 


The dynamic ground water model was run to predict the monthly 
variations in the heads throughout the area. This model used the 
transmissivity and mean recharge values determined in the steady state 
model and used the head distribution found from the steady state model 
as the initial heads. The dynamic model was used to estimate the storage 
coefficients throughout the basin, and to verify the variations with 
time of the recharge rates. 

The dynamic model was run for six years from October 1966 until 
September 1972, using a one month time step, and the calculated recharge 
rates (Figure 13). The model was also run for an additional twelve time 
steps using the long term mean infiltration rate to allow the heads to 
settle to the mean values before the real data were inserted. 

During initial runs of the model, using the small storage coef- 
ficient representing confined conditions, the modelled head variations 
at all the observation wells were considerably larger than those measured. 
The storage coefficients had to be increased to moderate the modelled 
head variations. The modelled head variations at the six observation 
wells were brought close to the measured heads by changing the storage 
coefficient around each well. Due to the small number of observation 
wells in the basin, the storage coefficients had to be changed over 
large areas around each well (Figure 33). The higher storage coefficients 
(from 0.04 to 0.16) used are consistent with water table conditions, and 
are Similar to the values derived by consideration of the ground water 
balance in Section 3.3.1. The storage coefficients derived from the 
pumping tests were not typical of the conditions over the majority of 
the basin. 

The modelled heads, using these modified storage coefficient 
figures, are compared with the measured heads in the six observation 
wells in figures 34 to 39. As the true elevations of the observation 
wells were not exactly known, and as a positional error was introduced 
when the wells were assigned to a node position, the elevations of the 
average modelled and average measured heads at the observation wells 
were somewhat different (Table 11). However, the variations in the 
modelled and measured heads generally agree fairly well at the six 
observation wells (figures 34-39). The standard deviations of the 
modelled and measured heads, listed in Table 11, agree closely at the 
Six wells, and the mean error between the head variations varies from 
0.53 feet at well 2B to 1.64 feet at well 4 (Table 11). 

The measured and modelled baseflows are compared in Figure 40. The 
baseflows generated by the model are the deep or slow responding ground 
water discharges (shown diagrammatically in Figure 21) derived from the 
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ground water reservoir extending over the whole basin. The measured 
baseflows, however, include the same deep or slow responding discharges 
as well as the shallow or fast responding discharges, some of which may 
be derived from local perched ground water reservoirs. In addition, a 
certain amount of interflow may be included in the measured baseflow 
figures, due to the method of analysis (Section 3.4.1). For this reason, 
the modelled baseflows are less variable than the measured values as 
shown in Figure 40. 


CONCLUS IONS 


The steady state model was found to simulate the measured heads in 
the basin closely with the average error between the measured and modelled 
heads being 6.5 feet and the largest error 28 feet. The total relief of 
the piezometric surface in the basin was 260 feet. In addition, it was 
postulated from the steady state model that the regional transmissivi- 
ties were about 1.5 times the mean transmissivities found from the water 
wells, and that the infiltration varied from 7 inches a year over the 
areas covered by more than 50 feet of till, to 15 inches a year over the 
areas near the creek where the till was thin or absent. 

The dynamic model simulated the head variations in the observation 
wells reasonably well. The storage coefficients around the wells were 
estimated by the requirements of the model response, and ranged from 
0.04 to 0.16 in different areas of the basin. The model simulated the 
deep seated baseflows, which were considerably less variable than the 
estimated baseflows. 
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Histogram of the differences between 
the measured and modelled mean heads 
for the nodes in the basin. 
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Figure 32. Comparison of the measured and modelled steady state heads in the 
Blue Springs Creek basin. 


Ee: 


tii tin =a Bie 


Ze 


IW 


Ss 


SS 
LX WON Ns IN 


NM NW ee nae 


N 

Sa 
eee ass 
. LY Z a 


Whine Pe 
ro 


pos 


Rae LLLa LALA 


Figure 33. Coefficient of storage values used in the model. 
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Figure 34. Modelled and measured heads at observation well BS -IA. 
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Figure 35. Modelled and measured heads at observation well 
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Measured water level in feet below ground surface 


Figure 36. Modelled and measured heads at observation well 
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Figure 37. Modelled and measured heads at observation well 
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Figure 38. Modelled and measured heads at observation well 
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Figure 39. Modelled and measured heads at observation well BS -5. 
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Figure 40. Modelled and measured baseflow values for the federal 
gauge sub — basin. 
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HY DROCHEMISTRY 


Chemical composition of water in a river basin is an important 
consideration in any hydrological study. The suitability of the ground 
or surface waters for use in industry, irrigation or for drinking 
purposes can be assessed, and ground water flow patterns can be esti- 
mated or confirmed by a study of the hydrochemistry. In karst basins 
where the landforms are largely developed by solution of the local rock 
strata, an analysis of the hydrochemistry is vital to an understanding 
of the karst processes. In particular, the relative saturations of the 
various minerals in the water are needed in order to assess the areas 
where solution of the rock is likely to take place. In this study, the 
hydrochemistry of the water was used to confirm the generalized ground 
water flow patterns and to study the saturations of the water with 
regard to dolomite and calcite in order to examine the karst development. 


WATER SAMPLING 


To investigate the lateral variations of the water chemistry over 
the Blue Springs Creek basin, 73 water samples were collected throughout 
the basin in July, 1971. Sixty-two samples were collected from existing 
bedrock wells and 11 samples were taken from springs and streams in and 
near the basin. The sample locations are shown in Figure 41. As there 
is only a small number of overburden wells in the basin (see Map 3), 
detailed sampling of the overburden aquifer was not possible. It is 
assumed, however, that the water samples taken from the overburden 
springs in the area, represent hydrochemical conditions in the over- 
burden aquifer itself. 

Samples of the well water were obtained from private homes using 
existing pumps and plumbing facilities to retrieve the samples from the 
wells. The water was allowed to flow for some time before a sample was 
taken to allow the stagnant water samples were stored in sealed glass 
bottles and transferred to the laboratory of the Ministry of the Environment, 
where the chemical analyses were carried out. The ions analyzed included: 
calcium, magnesium, sodium, potassium, chloride, sulphate, bicarbonate, 
nitrate and iron; pH and conductivity were also measured. The ionic 
concentrations expressed as milli-moles/litre (Hem, 1970) for the 73 
water samples are given in Appendix 2. 


POTABILITY OF THE WATER 


Drinking water should not contain constituents in large enough 
concentrations to render it harmful, unpleasant to drink or unsuitable 
for household purposes. Limits on the concentrations of constituents 
have been listed by the Ministry of the Environment (1974). The pub- 
lished criteria for private water supply are divided into two categories, 
the permissible and the desirable criteria. Water containing constitu- 
ents which generally exceed the permissible limits is not fully accept- 
able for water supply. Water that meets the desirable limits, which are 
more stringent than the permissible limits, provides good quality water 
Supplies. The permissible and desirable criteria for the constituents 
that were analyzed for in this study are shown in Table 12. 

tt was found that many of the water samples did not meet the 
permissible limits for total dissolved solids or for iron. Their locations 
are shown in Figure 42. Out of a total of 73 samples, 34 had total 
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with the Ministry of the Environment. These 
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and in these cases the depths for the wells 
were obtained by questioning well Owners 
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Figure 41. Location of water samples in the Blue Springs Creek basin. 
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TABLE 12. RECOMMENDED MAXIMUM PERMISSIBLE AND DESIRABLE 
CONCENTRATIONS OF SOME CONSTITUENTS FOR POTABLE WATER 
SUPPLIES, AFTER THE MINISTRY OF ENVIRONMENT (1974) 


Permissible criteria Desirable criteria 
(mg/1) (milli-moles/1) (mg/1) (milli-moles/1) 

Chloride 250 Tes <2) SOF L 
livatey a Ons - Virtuabery, absent 
Nitrate (as N) LO Oe | virtually absent 
Sulphate Z50 PREXS) =50 SOMo2 
Total Dissolved 
Solids 500 = <200 = 
pH 6..0-8. 5. unt cs 


TABLE 13. Classification of Hydrochemical Facies, after Back (1960) 


Percentage of constituents, in 
equivalents per million 


Ca+Mg NatK HCO 3 Cleat ToC, et NO3 


Cation facies: 
Calcium-magnesium..... 90-100 ag! AO rir Merk eS a iyo ae 
Calcoum-SoOGlume. oo ee: 50-90 LO=50% pow, co cceke cea sda otal eee eae 


Chloride-sulfate- 
lou Wer- Nalole 9% qos poy A Aue a IN ee Lae 10-50 50-90 
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Figure 42. Chemical potability of water samples in the Blue Springs 
Creek basin. 


83 


dissolved solids above the limit of 500 mg/l. Thirteen samples had an 
iron concentration above the 0.3 mg/l limit, although it is possible 
that the high iron content in some samples could be from the pumps and 
plumbing systems. In general, however, water can be treated readily in 
home systems to remove excess iron. 

Over half of the waters sampled in the Blue Springs Creek basin had 
one or more constituents which were above the permissible limits. 


Hardness of the Water 


Hardness is a measure of the amount of calcium and magnesium in 
solution. Hard water is objectionable because it tends to form scale in 
water heaters, pipes, kettles and boilers, and because a large amount of 
soap is needed to form a lather. Water hardness is often classified from 
soft (less than 60 mg/l equivalent CaCO3) to very hard (over 180 mg/1 
equivalent CaCO3), (Dunfor and Becker, 1964). All the water samples in 
the Blue Springs Creek basin (except sample 60, which is discussed in 
Section 5.3) were very hard (Appendix 2), containing from 216 to 600 
mg/l of equivalent CaCO3. The majority of this hardness is of the 
temporary type (e.g. carbonate hardness) and could be reduced by suitable 
home treatment systems. 


GENERAL CLASSIFICATION OF THE WATER 


In order to assess the differences in the chemistry of the surface 
and ground waters in the basin, the water samples were classified into 
hydrochemical facies using the terminology developed by Back (1960). 
This classification is based on the dominant anions and cations present 
in the water. The various hydrochemical facies are listed in Table 13 
and have been inserted on the hydrochemical facies diagram, Figure 43. 
Generally, a natural water will range from a predominantly bicarbonate 
water, through a sulphate water, to a predominantly chloride water as 
the degree of 'metamorphism' of the water increases (Chebotarev, 1955). 

The percentages of milli-equivalents per litre (meq/l) of the 
various cations and anions were plotted on the trilinear and parallelo- 
gram diagrams shown in Figure 43. It can be seen from Figure 43 that 
the majority of the water samples were calcium-magnesium, and bicarbonate- 
chloride-sulphate facies. Most of the water is therefore in its first 
or second phase of 'metamorphism'. 

Sample 42 is unusual in that the potassium concentration is very 
high. Potassium rarely exceeds the sodium concentration (Hem, 1970), 
but in sample 42 the potassium concentration is over five times the 
sodium concentration. The chloride, sulphate and nitrate concentrations 


are also high, leading one to suspect that this well is being contaminated 
by fertilizers, sewage or farmyard wastes. 


Areal Variations in the Degree of Metamorphism of the Water 


The lateral variations in the degree of metamorphism were investi- 
gated. The chloride and sulphate to bicarbonate ratios were plotted far 
the basin as shown in Figure 44 (in pocket), where the high numbers 
indicate a relatively high degree of metamorphism. 

There are some general trends in Figure 44, the predominant one 
being the banding roughly parallel to the creek in the middle and 
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Figure 43. Chemical composition of water in the Blue Springs Creek basin. 
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southerly parts of the basin. The degree of metamorphism appears low in 
the ground water recharge and discharge areas and higher in the areas 
between. In a general way the degree of metamorphism over the basin 
fits the ground water divide (on Map 7) more closely than the surface 
water divide. 


DEGREE OF SATURATION OF WATER WITH RESPECT TO DOLOMITE AND CALCITE 


Chemical solution of the aquifer material will play an important 
part in the erosion and geomorphic development of any areas underlain by 
carbonate rocks. In karst areas in particular, the chemical erosion of 
the rocks predominates over mechanical erosion. The solution of the 
rocks by the water is limited by the solubility of the aquifer material 
in the water. Thus an understanding of the solubility of the aquifer 
material is important for the understanding of the geomorphic develop- 
ment of carbonate areas, and to assess regions where solution of the 
aquifer material is likely to be significant. 

Water which is in contact with soils and carbonate rocks will 
dissolve some of the minerals. Given time, the solution of the minerals 
will continue until the water becomes saturated with respect to that 
mineral. When this occurs, the solubility equation is applicable; for 
dolomite this can be expressed as: 


bh ++ -- 
iCal ue Mg wae COs me). 
K Hs ETO) 
{Ca Mg (CO3) 2} 
++, Seer . : : 
where {ca } is the activity of the calcium in solution; 
++, aa : i 
{Mg } is the activity of the magnesium in solution; 
1COnmmn: is the activity of the carbonate in solution; 
{ca Mg (CO3) 24} is the activity of the solid dolomite (which 
by convention is unity); 
and Ky is the solubility product for dolomite (which 
is temperature and pressure dependent) . 
Equation 17 only applies if the water is in contact and in equilibrium 
with solid dolomite. 
If the water is not saturated with respect to dolomite, the water 
1S considered to be aggressive and is capable of dissolving dolomite. 
Conversely, the water can be Supersaturated, particularly if carbon 


dioxide is lost, and precipitation of dolomite will tend to occur. A 


convenient index for the degree of saturation can be given by (Wigley, 
19/1)" 


SATDOL = log ({ca’*} {g' } {COs i) ellog {x} (18) 
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where SATDOL is the saturation index with respect to dolomite: it is 
positive if the water is supersaturated with respect to dolomite and 
negative if undersaturated. Comparing equation 17 with Equation 18, it 
can be seen that SATDOL will be zero if the water is saturated. 

A similar index can be given with respect to calcite as: 


PATeAL = Log (ca } 4c0,.° 1 — log (Kk) (19) 


where SATCAL is the saturation index with respect to calcite, 
and K, wertcne sOlubijaty product for calcite. 

In order to determine the saturation indices, the activities of the 
appropriate ions are needed. The activities and saturation indices can 
be derived from the chemical analysis figures. The calculations are 
more complicated if ion pairs are considered, (Garrels and Christ, 

1965). Ion pairs have often been neglected in hydrochemistry calculations 
(for example Van Everdingen 1969), but Wigley (1971) has shown that ion 
pairs can be significant in hydrochemistry calculations and appreciably 
alter the calculated ion concentrations and degrees of saturation. In 
this study, the following ion pairs were considered: CaS0,°, MgSO,°, 
CaCO,” , MgHCO3*, Cac03°, \MgCO39, NaSO,7, KSO, , NaCO3-, NaHCO3°. 

Other possible ion pairs (for example KCO3 , CaCl*t) were considered not 

to be present in significant amounts. 

The calculation of the ion and ion pair concentrations was carried 
out using a computer program which was written and has been described by 
Wigley (1971). The program calculates the degree of saturation with respect 
to calcite and dolomite; the results are listed in Appendix 2. 


Errors that Affect the Calculated Degree of Saturation of Water with 
Respect to Calcite and Dolomite 


There are a number of errors, which are associated with the sampling 
and the analysis of water, that may affect the calculated degree of 
saturation with respect to carbonate minerals. These errors and their 
effect on the calculated degree of saturation are described in this 
section. 

Certain chemical changes may take place in the water between the 
Sampling and analysis time (Hem, 1970). If the water is supersaturated 
with respect to any mineral, precipitation will probably take place to 
bring the water towards equilibrium. Precipitation of calcite or dolomite 
will decrease the measured amount of calcium, magnesium and alkalinity 
and will also decrease the measured pH. 

If the water is undersaturated, it is possible to lose dissolved 
carbon dioxide from the water if the partial pressure of the carbon 
dioxide is above the atmospheric value, particularly if the sample 
bottles are not completely filled, or if the bottles are left open for a 
time before analysis. Loss of carbon dioxide will result in a decrease 
in the alkalinity and an increase in the pH value, and indirectly will 
result in a reduced degree of calculated undersaturation with respect to 
calcite or dolomite. If sufficient carbon dioxide is lost, the water 
may become supersaturated and start to deposit carbonates. To determine 
the changes in the water, the most sensitive parameter to be examined is 
the pH: if the pH rises, the water is losing carbon dioxide, and if the 
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pH falls, carbonates are probably being deposited without loss of carbon 
dioxide. In many natural waters, the pH rises during sample storage, 
indicating that carbon dioxide is being lost (Hem, 1970). 

To estimate these changes in chemistry that may occur in the water 
samples from Blue Springs Creek, three samples were collected on 
December 17, 1973, at the sites marked on Figure 41. The pH and temper- 
ature were measured in the field and subsequently analyses were performed 
in the laboratories of the Ministry of the Environment. The relevant 
data are shown in Table 14 and Appendix 2. It can be seen that in all 
three samples the pH did rise somewhat between sampling and laboratory 
analysis likely due to loss of carbon dioxide from the samples during 
storage. The calculated partial pressure of carbon dioxide is listed in 
Table 14 and can be seen to be well above the normal atmospheric partial 
pressure of 0.03 per cent. The largest pH change occurred in sample 83, 
which had the highest partial pressure of carbon dioxide and would 
therefore be most likely to lose the largest amount of carbon dioxide. 

The effect of the temperature changes can also be seen from Table 
14. The water sampled varied from 7.0 to 8.5°C and the laboratory was 
assumed to be at 20°C (68°F) during analysis. As water is warmed, 
carbonates become more soluble and the degree of saturation is reduced. 
Samples 81 and 82 (which were from surface streams) were both supersatur- 
ated with respect to calcite at the time of sampling. When analysed in 
the laboratory, all samples had equilibriated and were all approximately 
Saturated with respect to calcite. The degrees of saturation with 
respect to dolomite were generally lower, although samples 81 and 82 
were supersaturated at the sampling time. 

Generally speaking, the saturations calculated using field pH 
Measurements are below those using laboratory pH measurements. Thus the 
Saturations in Appendix B, which use the laboratory pH and an assumed 
field water temperature of 10°C, are probably too high. Changes in the 
degree of calcite and dolomite saturations should, however, be generally 
less than 0.3 and 0.6 units, respectively. 

Under ideal conditions, pH measurements can be accurate to 0.05 
units. Thus, as the calcite saturation is approximately proportional to 
the hydrogen ion concentration, and the dolomite saturation approximately 
proportional to the square of the hydrogen ion concentration, the corres- 
ponding errors in the degree of calcite and dolomite saturations are 
about 0.05 and 0.1 units, respectively. Thus, an analysis which had a 
degree Of calcite saturation trom —0.05 to) 0.30) units, sora egreesor 


dolomite saturation from -0.1 to 0.60 units, was considered to be fully 
saturated. 


Areal Saturations 


The degree of saturation of the water determines if the water is 
capable of dissolving the minerals in the aquifer material. If the 
water is undersaturated, solution can occur. If Saturated, the water is 
in equilibrium with the minerals and no solution or precipitation of the 
minerals will occur. If the water is supersaturated, which can occur 
when ground water resurges and loses carbon dioxide to the air, preci- 
pitation of minerals is probable. 

The degree of saturation of the water with respect to calcite and 
dolomite was studied in the Blue Springs Creek basin to determine the 
areas where solution and precipitation are probable. The degree of 
Saturation of the water samples is shown in Figure 45 (in pocket). Most 
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of the samples (47 out of the total of 73) were saturated with respect 
to calcite. Only four samples were supersaturated, and 22 samples were 
undersaturated. Most of the samples (55) were undersaturated with 
respect to dolomite. Only three samples were supersaturated and 15 were 
saturated with respect to dolomite, although all the well samples were 
obtained from dolomitic bedrock wells. The reason for the lack of 
equilibrium with respect to dolomite is further discussed in Section 
Seo. Os 

In Figure 45, it can be seen that there are no obvious trends in 
the degree of saturation with respect to dolomite. The undersaturated 
samples appear to be randomly distributed. The reason for this effect 
teoeexplained in Section, 5.5.3. 

The calcite figures also appear to be fairly random; however, many 
of the undersaturated samples are near the drainage divide where the 
water is recharging the aquifer. Conversely, many of the saturated 
samples are in the centre of the basin, where the ground water has had a 
long residence time in the aquifer and so has had time to equilibriate 
with the aquifer material. 

When the ground water resurges it generally gives off carbon dioxide, 
as the partial pressure of carbon dioxide is usually above that of the 
atmosphere. This increases the degree of saturation of the water and 
may give water which is supersaturated with respect to calcite and 
dolomite. In this case the minerals would tend to precipitate out. If 
the ground water resurges at a bedrock spring the precipitate of the 
carbonate can lead to a deposit of tufa; such deposits were seen at 
Rockwood just to the west of the basin. In the Blue Springs Creek 
basin, however, the only water which is supersaturated with respect to 
both calcite and dolomite is in the main stream (Figure 45). Any preci- 
pitation that occurs would be in the sediments giving a calcareous mud 
or marl deposit. 


CALCIUM TO MAGNESIUM RATIO 


The calcium to magnesium ratios (Ca/Mg) of the waters in the Blue 
Springs Creek basin were used to estimate generalized ground water flow 
patterns in the bedrock aquifer. In general, ground water will tend to 
equilibriate with the surrounding rock, and with time, the (Ca/Mg) of 
the water will approach the (Ca/Mg) of the aquifer material. As the 
chemical changes are not instantaneous, the water which has had a short 
residence time in the bedrock aquifer, will have a different (Ca/Mg) 
than the water which has had a long residence time in the bedrock aquifer. 
Generally the ground water which is being recharged has had a shorter 
residence time than the ground water that is being discharged. It was 
found possible to predict, in a general way, the recharge and discharge 
areas of the bedrock aquifer by using the calcium to magnesium ratio of 
the waters. To do this it was first necessary to determine the (Ca/Mg) 


in the bedrock and overburden materials; this is described in the fol- 
lowing subsections. 


Calcium to Magnesium Ratio of the Amabel Dolomite 
The (Ca/Mg) of the Amabel bedrock is given by rock analysis figures 


quoted by Hewitt (1964). A mean of six samples from quarries near the 
towns of Milton and Guelph (the quarry localities are shown on Figure 2) 
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gave a CaO and MgO composition of 30.71 per cent and 20.19 per cent, 
respectively by weight. Converted to the mole fraction (a mole is the 
molecular weight of substance in grams), the calcium/magnesium mole 
Tatwo 1S? 


(Ca) eee 087A 40.31 
(Mg) 56.08 20.19 


r.09 


as the molecular weights of Ca) and Mg) are 56.08 and 40.31, respectively. 
The extreme values for the seven samples gave a (Ca/Mg) ratio of 1.01 at 
Guelph yto7l.730 atl Militon.,, Cowell and Gregor (1973) found a molar ratio 
from 1.00 to 1.08 with a mean of 1.03 in five samples of the Colpoy Bay 
and Wiarton members of the Amabel at Rockwood. 

The (Ca/Mg) of pure dolomite (CaMg(CO3)9) is 1.0. Thus the Amabel 
is a fairly pure dolomite and contains only a small percentage of 
calcite. 


Calcium to Magnesium Ratio of the Wentworth Till 


The carbonate composition of the Wentworth till in the Guelph area 
is quoted by Karrow (1968). The calcite/dolomite ratio was measured 
uSing the method of Dremanis (1962), and the average ratio of 25 samples 
of Wentworth till was 0.63 by weight. Converted to the mole fraction, 
Ene ratio Ls: 


(Calcite) 07.63 184.41 

== =S = xy = = 1.15 

pe torcyee ear T00n0o? “4 Tal 

as the molecular weights of calcite and dolomite are 100.09 and 184.41, 
respectively. As one mole of dolomite contains one mole of calcium and 


one mole of magnesium, the (Ca/Mg) of the till will be: 


(Ca) Is i5+1,0 

ae Se = 2.15 

(Mg) mole tae 

The extreme values for the (Ca/Mg) of the 25 samples of Wentworth 


Prieare. ie 5o-and 2.96. 


Comparison of the Calcium to Magnesium Ratios of the Water and of the 


Aquifer Material 


The molar ratios of the dolomite bedrock and overburden have been 
added to the trilinear diagram, Figure 43. The majority of the waters 
have a (Ca/Mg) similar to that of the overburden material, although all 
the wells were in the dolomite bedrock. 

In Pennsylvania, Miesler and Becker (1967) found a close correlation 
between the (Ca/Mg) of the ground water and the bedrock lithology. 

Wells in dolomite gave a mean (Ca/Mg) of 1.09, while limestone wells 
gave a mean ratio of about 4. However, the area in Pennsylvania was Loe 
covered with much overburden, while in Blue Springs Creek the bedrock is 
covered with a considerable thickness of till. The till must be affect- 
ing the (Ca/Mg) in the underlying dolomitic aquifer. 

In the Blue Springs Creek area the water in the till generally 
becomes almost saturated with respect to calcite (as can be seen from 
samples 39 and 7, Appendix 2, which are from springs in the CL LE ee 
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is undersaturated with respect to the dolomite. When this water enters 
the dolomite bedrock, solution takes place, but as the dolomite dissolves, 
the calcite concentration must increase and the water becomes supersatur- 
ated and thus precipitates calcite. This calcite deposition hinders 
dolomite solution and prevents the water from reaching saturation with 
respect to the dolomite and also from approaching the (Ca/Mg) of the 
dolomite bedrock. 


Areal Variations of the Calcium to Magnesium Ratio of the Water 


The (Ca/Mg) of the waters throughout the Blue Springs Creek basin 
were plotted in Figure 46 (in pocket). The ratios appear higher near 
the basin divides and become lower towards the centre of the basin. The 
areas of high (Ca/Mg) represent areas of bedrock recharge and the areas 
of lower ratio are areas where the water has had a long residence time 
in the dolomite; they also tend to be areas where the bedrock is dis- 
charging water. In the extreme northern part of the basin, the (Ca/Mg) 
is low, which supports the findings of Section 3.3 that this area is a 
discharge zone for the bedrock aquifer, in spite of the fact that the 
area is near the topographic divide. 
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SUMMARY AND CONCLUSIONS 


A hydrogeologic study was undertaken in the Blue Springs Creek 
basin to evaluate the ground water conditions in terms of distribution, 
quantity and quality. The basin is representative of karst conditions 
in Southern Ontario, such as those found along the top of the Niagara 
Escarpment. A summary of the significant results of the study, and the 
conclusions derived from these results are given below. 


1) 


2) 


3) 


4) 


The hydrogeologic studies reveal that the materials underlying 
the Blue Springs Creek basin form a complex hydrologic system. 
This system comprises a dolomite unit, a sandy till unit anda 
kame, Outwash and alluvial deposit unit. The dolomite unit 
(Amabel dolomite) forms the principal aquifer in the study 
area. The kame, outwash and alluvial deposit unit combined 
with the sandy till forms a minor aquifer. 

Most wells constructed in the dolomite bedrock aquifer 
should yield adequate quantities of water for domestic farm, 
commercial and smaller industrial uses, but the water will 
likely be very hard and have a total dissolved solids content 
above the 500 mg/l permissible limit recommended for public 
and private water supplies. 

Based on geomorphic and hydrogeologic evidence, the develop- 
ment of karst features in the Amabel dolomite is minor. There 
are several caves and springs at Rockwood and a number of 
surface sinks and springs are found in the basin which interact 
with the dolomite aquifer. As the karst is not well developed, 
it does not play a major role in altering the ground water 
flow patterns. 

Conditions during the study period (1967-1971) were found to 
be slightly wetter than the long term mean conditions. The 
precipitation at Guelph during the study period averaged 33.7 
inches/yr (855 mm/yr), which was 3 per cent higher than the 
ToAI=VO7Osmeantor 32.7 inches/yr(S30°mm/yr). “The streamflow 
on the Speed River below Guelph (at the federal gauge O2GAO15) 
averaged 212 cfs (6.0 m3/s) during the study period, which is 
8 per cent higher than the mean flow for the period of record 
which was 197 cfs (5.6 m3/s). 

The transmissivity of the overburden was estimated from well 
records and the median transmissivity was found to be 3,500 
IGPD/ft. (600 mm?/s). A large range of transmissivity values 
was found; the 10 and 90 percentile figures being 300 and 
42,000 IGPD/£t (50 and 7200 mm?/s). 

The transmissivity of the bedrock aquifer was estimated 
using three different approaches. Calculations were carried 
out using pumping test data, well record data and calculated 
ground water flow rates. The average transmissivity values 
found by the three methods were 3,500, 3,200 and 6,950 IGPD/ft 
(600, 550 and 1200 mm2/s) respectively. The 10 and 90 per- 
centile figures for the bedrock transmissivity determined from 
well record data were 610 and 17,000 IGPD/ft (105 and 3000 
mm“/s), a significantly smaller range than was found for the 
overburden wells. 

No significant differences in the transmissivity values 
were found when the bedrock wells were divided into classes 
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5) 


7) 


8) 


2) 


10) 


according to the static water levels. It was found, however, 
that the wells where the aquifer was shallow showed a signi- 
ficantly larger range in transmissivity than the deeper wells. 
This supports the idea of shallow development of karst solution 
channels and conduits. 
The extent of the ground water basin was estimated from the 
piezometric map and was found to have an area of 26.4 square 
miles (68.5 km“) which is 13 per cent smaller than the topo- 
graphic basin area (30.2 square miles or 77.3 km?) . The small 
difference between the topographic and ground water basin 
areas supports the theory that the karst is not well developed. 
The bedrock aquifer was found to be recharged at about 10 
inches (254 mm) per year, primarily by the overburden aquifer. 
Discharge from the bedrock aquifer is at discrete springs 
along the creeks in a number of swamps, and along portions of 
the creek beds. 
The surface runoff component of the streamflow was separated 
from the ground water discharge component using conventional 
baseflow separation techniques. The baseflow was found to 
average 14.0 cfs (0.40 m3/s) during 1966-69 at the federal 
gauge, which is approximately 68 per cent of the streamflow of 
20.6 cfs (0.58 m3/s). The average baseflow was equivalent to 
an average annual ground water runoff rate of 9.8 inches/yr 
(250 mm/yr) from the area of the basin. 

The relatively high degree of consistency of the reces- 
Sion curve determined from the stream hydrographs, indicates 
that the majority of the ground water discharge comes from one 
principal aquifer or several aquifers having similar hydraulic 
characteristics. 
The evaluation of the specific yield in the different sub- 
basins during periods when recharge by precipitation and 
evapotranspiration are assumed negligible, can be used as a 
measure of the characteristics of the aquifers feeding the 
streams. The yields vary from 7.2 per cent in the Boy Scouts 
Camp sub-basin, to 12.7 per cent in the Eden Mills sub-basin, 
indicating that the yields in the upstream part of the Blue 
Springs Creek basin are smaller than in the lower sections. 
The higher yields were related to large areas of thick persis- 
tent kame, esker and outwash sands and gravels, which occur 
principally in the lower sections of the basin. 
The steady state piezometric surface was successfully modelled 
in the basin using finite difference techniques. In a total 
piezometric relief of 280 feet (85 metres), the heads were 
Simulated to an average error between the modelled and measured 
heads of 6.8 feet (2.1 metres) with the largest error being 28 
feet (8.5 metres). The transmissivities used in the model 
were 1.5 times larger than the median transmissivities found 
by well analysis. 
Incorporating variations with time, the ground water flow 
regime was modelled on a monthly basis from 1966 to 1972. The 
variations in the heads at the six observation wells in the 
basin were modelled accurately, with the mean difference 
between the modelled and measured head variations being 0.95 
feet (0.3 metres). The storage coefficients throughout the 


basin ranged from 0.04 to 0.16. The recharge to the ground 
water, determined by water budget considerations, was verified 
by the model. 

11) Expe~ience with the ground water models showed that detailed 
input data are required to simulate the piezometric heads 
satisfactorily. It was found that the results of the models 
are sensitive to changes in transmissivity, coefficient of 
storage and ground water recharge values used in the model. 

It was difficult to generalize the parameters in the Blue 
Springs Creek representative basin, as all of the sensitive 
parameters were found to vary considerably across the basin; 
the ground water recharge also varies significantly with time. 
The only data not sensitive to the model operation are the 
Yiver permeabilities; high values can generally be used for 
the modelling. Thus, in examining the the feasibility of 
extrapolating the Blue Springs Creek model to other areas of 
Similar physiography, basic data on the three sensitive input 
parameters would be required for the area to be modelled. 

12) Thirty-four of a total of 73 ground water and surface water 
samples were found to have a total dissolved solids concentra- 
tion above the permissible limit (500 mg/l) recommended for 
drinking water. All except one of the samples were found to 
be very hard, although much of the hardness appears temporary 
and could be reduced easily by suitable home treatment facilities. 

13) The majority of the water was found to be of the calcium- 
magnesium and bicarbonate-chloride-sulphate facies type. Two 
samples were found to be of different hydrochemical facies; 
One is likely due to cation exchange through the clays and the 
other possibly due to contamination by fertilizers or sewage. 

The relative degree of metamorphism of the water was 
found to vary across the basin with the lowest metamorphism 
Occurring near the recharge and discharge zones. Most of the 
water samples are in the first or second stage of metamorphism. 

14) The ground water was generally found to be saturated with 
respect to calcite and undersaturated with respect to dolomite. 
Surface waters were often supersaturated with respect to both 
minerals. The surface water would therefore tend to preci- 
pitate calcite and dolomite. 

15) The calcium to magnesium ratios (Ca/Mg) of the waters were 
compared to the (Ca/Mg) of the aquifer materials. It was 
found that the waters had ratios comparable to the overburden 
aquifer material, even the samples from the bedrock wells. 

The water infiltrating down through the overburden becomes 
saturated with calcite and when the water enters the bedrock 
it is unable to equilibriate with the dolomite, due to the 
deposition of calcite in the aquifer material. Thus the water 
never becomes saturated with respect to dolomite and never 
approaches the (Ca/Mg) of the dolomite. 

The areal trends in the ratios were investigated. Flow 
patterns of the ground water can be inferred from the (Ca/Mg) 
trends; the flow pattern was found to match, in a general way, 
the pattern found by consideration of the piezometric surface. 
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16) 


The trends in water chemistry match the ground water flow 
patterns and ground water divide location better than they 
Match the topographic flow patterns and topographic basin 
divide, confirming the ground water flow patterns found by 
consideration of the piezometric surface. 
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APPENDIX 1 
Water Budget in the Blue Springs Creek Basin, 1966-1972 


A summary of this Appendix is given in Table 6. 
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APPENDIX 2 
Chemical Composition of the Water Samples 


in the Blue Springs Creek Basin 


VLG dina § Gr Wey OG0ORS Oe) = OSE 00°0 It S@ GORO SO Ov*O Mb iy COME Syhe 1h SE 


GIES GOm iS ES Hid O = Oi ORS OSL SL 70) Oesr ELS O S\~O) ES @ v8°s Scar SV aEG vE 
OVE COs 9°60S Oc*O 56> O) OLE L0°O 8c°O €0°O EieO i? 0) GES LOaT [SEG ce 
897 SOs ELS: MLO = LO SOnS OSs~z S SiO) 6L°O SOG VLaO 66°0O CS) Semel SSec GE 
oT Com wwe Cie Om 80°0 OS EO) Sian) CEC 9Gs0 CS) Ik ILS SG 1 GOmG 1S 
STE cO* T° S60 hs 0) = SOsORS OSm2 GEO ere SGmO €8°0 Ovo 80°0 SHE UE cO°? O€ 
OL OT* 9°8CU O€*O O€*O OomZ LORO cv Oo €0°O O€*O OSe0 Comvy, c8°O S8°T 6¢ 
80E SO° Sai? ss) = LOO OSes SOMO UE oO €0°O IE~@ IKE @ OLS 98°0 OGaG 8 
OVE SL iy BES SO = vo°o OSE GO Es) FC SEsO Lore v9°O XS c8°O GSEG LG 
86C Om 9° 0D Ov O = SO" = OST 60°0 0c °O €0°O Se SEO GL Ws 98°0 Glues 9g 
OTE Oy SROLY, Sit) = 80°0 O9°L 9050 les, @) S0)~O) IEE cv O 88°77 S6°0 CUrG SiG 
Lh SO* 6° 789 SiO) < 6GORORS Oiler 00°O cv 0 Lio) O€*O he) CSE OV'T Gmc vc 
6LC Site OEY SO = 60°0O OSE 00°0O SiGe O SORO nL) c0°O SIe~ c8°O Soir Ec 
86C ie Si GSiz SiO = TO°O 0992 00°0 82°O €0°O 6G O€*O 89°70 Ue ESaat KG 
8LC cals 9° TDP QQ = €0°O O9g2 00°O SGaO S0*O ClO 82°O 00°S 66°0 08 °T LG 
EVE Oy DQ ES iz Oy = LO-O = Ores SOmO is SO0°*O O€*O 6S°0O OGES Dice ASG 0c 
REE SO)” 8°8TS EOmO mS iLO OE SORO CIO 80°0 Sean0 cv-O Gems EKO) 1 OeeG 6T 
VC O° 8° 68E NG) = €0°O OLE 00°0 Sv °O €0°O OO) 6c°O EM VEAO OLSE 8T 
STE gir” c°Tev Si-O = OE lS OV°L S0°O 87 °O SO~ is 7) vv-O VO°S 66°0 OGG a 
TSv SO” UL So @ = SOO = OG LE’ L8°0O OT*O T6°0O 69°0O 80°9 Slvr CASAS OE 
GES SO* 8°¢80 EO = LOMO! = OSRL 90°0 Ovo vo°O GemO 60°O 88° S6°0 fet @ Sip 
VE Sv* ORES SOnORS SiO O9Rz LE @ c8°O vo°o 6€°O Ovo 86°7 S6°0 OS.aG vT 
Gas SO° f° 680 SO = SEO S Oe ScanO Sao 60°O 1 -O) Sica O 88°70 T6°0 CIRC or 
OVE SOx iv ES O)"@ = LO) = Ocee Sie O Bea) vI“O O€*O Gir Oc*S We AE SoG Gil 
abe 1 L°8Lv OG me c0°O OS. cO°*O a ae) SOMO 60°0 Nee CSS, EOmE SOT 1248 
SRELS SO° Si GES; Viva On LOSORS Ov~L 00°O EGO c0°O EO S950 WES SO. AS~G OT 
EC Sie © LGV, 82°O como) 06°L TO°O 87 °O €0°O CEO LGW 09°? vL°O Omit 60 
BLS SS)” OSs GO = KOO) = OLE 00°0 90°0O €0°O OSE EO) WG 98°0 OfSL 80 
Ge SO 6 SiLiy SemOm LOO OSael: cT*O cv O TO0°O alee O Tv°O 88°0 c8°0O OfnG LO 
COV ibe ie v3 QL" O = Or = OCae ke COM €0°O vO°Tt 08°0O 78°S c8°O O8°e 90 
VCS SO 1 EAS} AZO) c0°O OGere Ov ~O GE IE 7O°O 00°? LOO VIS Cileal SO°? SO 
OVE Tes SSIES [Ei Or = SOM. Ov*L SO°O Sas0 0c °O EiaO 08 °O vO°S cOnE INS vO 
Lev SO* EGE OE "O@ = do @ = Osa e EGO) 89°O ce °O Som SL 8) NES OVE LSZG €0 
80€ SO He Whshiiy <7 0) = CieOn = OV*L €0°O Ov*O OT*O O€*O ie OGsS eye ae EOmL cO 
€O€ ql S°8S0 IEE @) vo°o OO L TO0°O OO 70°O SIL EvmO 96° Locate GOnE TO 


oe) T/bur T/bu 9s 3TwoTod 83TOTeOD H SON nO) Si WN DEHOS _£OOH Bw ize) 
= = + an ++ ++ ord 
T/bu se uorzlL SPTTOS Oo} get ON 9eTaues 
ssoupizeH PeATOSSTC qZoedsezx yzTM uot A SAVTT/SETOW-TTTtTwW UT UOTRJeEAIUSOUOD TeROL 
TeqIOL -eznjies Jo s0769q 


Se ee ag aes oe ee eee ee 


106 


(0/0) 18) 
T/bu se 
ssoupieHyH 


SO° 
60) 
SO° 
Sim 
EO” 
SOs 
t- 
Gow 
SOs 
Om 
9° 
Sia 
6° 
SL 
SOs 
SOM 
Om 
SOs 
T° 
‘SKO)" 
SO° 
SO° 
t- 
[- 
SOm 
T° 
ze 
t 
SOY 
SO° 
Oy 
SL 
Sis 
GOR 


T/bu 
uOIT 


Z°T6P 
Z°O€S 
8°S99 
L°T8S 
S°9Sb 
T°S@p 
T° OSD 
€°6LY 
T°18p 
€° 709 
O°L6E 
L°ESE 
€°6SE 
8°L8S 
Z*Z0S 
6° Lt 
S°8tP 
9°6E€S 
bes 
6°80S 
T° ver 
€°P8P 
L°8@p 
6°L9S 
T°vEes 
6°79S 
Z°L8P 
LivLVY 
L°0S9 
9°66S 
L°080 
CELE 
6°S8E 
6°€8S 
T/bu 
SPTTOS 
PeaAToOsstaq 
TeOL 


67°00 - €0°0 
vr'O - 2170 
cf°O - 90°0 
Or'oO - v0°O 
€Z7°0 ECO 
c£°O - Z0°O 
SOL €9°0 
Ov'o - S0°0 
ZT*O - one) 
LDCE - RT 
vZ°O - 00°0 
LEO - 60°0 
6T°O - 00°0 
€0°0 - €T°O 
To"O - C10 
€0°O - vT°O 
€0°0 ST“O 
LE 700 = 90°0 
eT O01 80°0O 
6T°O - 90°0 
co (olsaae) 
€0°0 21-0 
06°0 8S°0 
9S°0 - Ste) 
O€"0 - £O-0 
00°0 - LUO 
ZE°O - S0°0 
SE°0 ore 0) 
ZT°O - L0°0 
50 — Z0°O 
Se°O - 90°0 
€£°O - 80°0 
LEO - TL 20 
LE°O o¢ "0 
2} TWOTOC @3TOTeD 


o2 
qoedsez yqytmM uot 
-einjzes fo sver7beq 


OS°L 
ONE Ue 
OV L 
OV -L 
OSiee 
OS*L 
07°83 
OV'L 
OS*L 
OLRL 
09°L 
Ore 
OLeE 
OS=Z 
ONL 
OL -/ 
OLee 
OR ee 
OD 
OS*L 
OGee 
Xi) ae 
02°83 
Ole“, 
Ola 
09°L 
OS*L 
OV-L 
OV°L 
OV°L 
OS*L 
O9mE 
O09°L 
OLS 

Hd 
get 


qo 
‘S) 
je) 


On 
(aoe Soy Lean aN ee oat | 
je) 


. 


Ses 


SLOP OCOLOPOOFOs OT Ore@ 


Ors LO i i OO 
eH ee qo 


jo) 


SO°O CEO 72°0 00°S IGE CTE 
v0°O LEO veo 96°S S Com 
80°0O T9°T 80°0 Gila S6°0 
9S°0 LS°O 99°0 9E°S (KO) 
70°O SLO LEC BL°P S6°0 
70°O SilanO €0°0O 08°? S6°0 
€0°0 S70 ClO 98° WES 
€0°O GE°O 80°0 8Z°S SO 
OG °O O€°O L0°O 9T°S Gomir 
TO°O GL°S BS % 7S°S Z0°O 
COnO Liew vo°O 80° S6°0 
COE SiO €0°0 [f° 98°0 
z0°O oc °O LOZO Comat 98°0 
€0°O OL°O 0S°0 00°9 61° T 
Z0°O 97°00 0€°O O7°S WIE 
SO0°O Sat 67°O Cli 98°0 
z0°O eae) veo 96°97 S6°0 
OT°O CEO 69°0O TE°S 61 °T 
ESE vv°O €£°O 88° 66°0 
LE o@ Aten O TS°0 vps 9E°T 
SO°O ale © 67°0 8L°P LO? 
€0°0O vL°O Se°O 9L°? LOT 
€0°0O 97°0 OS °O p9°D 06°0 
SOO T9°0 8S°0 90°S iTalent 
vL°O S9°0 67 °0O 0Ss°8 S61 
80°0O OL°*O LS°0 86° we 
vZ°0O o£ °O 67 °O 80°S SOmar 
BE °9 CE OO°E€ €v'6 97°? 
Ty °0O 6£°0O zZ9°O ZL°S CEO 
€0°0O 00°T mo OL°S €O°T 
TO°O Er oC 87°0 0z2°S vL°O 
Z70°O 9Z°0 Z0°O 89°P7 98°0 
€0°0O 9Z°0 0€°O VED S6°0 
70°O 8P°O LE°C 00°9 87°T 

Bi UN BtoOS _©OOH W 
+ + ++ 


SAATT/SETOU-TTTTW UT UOTJeEATQUeDUOD TeIOL 


SO 

eC 

L6 

Sv 

SE 

OT 

c6 
EG 
SE 
cO 
c8 
Lv 
Lv 
IS 
SG 
S6 
c6 
Sis 
ST 
GE 
08 
c6 
08 
LV 
Oc 
Sil 
(as) 
GE 
All 
09 
LES 
cv 
Sv 
SU 


++ 


XG 
nc 
eG 
NG 
KC 
aC 
aL 
XG 
NG 
“O 
aU 
mas 
pa 
<G 
nc 
‘Al 
a 
“G 
G 
wale 
a 8 
romp 
a 
<G 
ee 
G 
a2 
aS 
Te 
ac 
eG 
“AG 
aE 
%G 


170) 


69 
89 
ES 
NS, 
S9 
v9 
re) 
(Be, 
We 
09 
6S 
8S 
ES 
MS 
SIS) 
vS 
ea 
éS 
"S. 
OS 
60 
80 
Lv 
90 
Sv 
vv 
ev 
cD 
Tv 
OV 
Og 
8E 
Le 
9E 


ON oe Tdures 


107 


"€p eanbtq uT uMmoys 
aze soTdues 7eQeM SYR JO SOANOS pue UOCTREOOT SUL “*P°S UOTIOSS UT peqTAOSEep ST eqaTWOTOP pUe sqTOTeO 03 4Oedsez YATM UOT eEANFeS FO serzbep suL 


(ee SE SS SS ee ee 


86C¢ SO° 6° OED Ov°O 9620 Gee Ew S3e0 €0°O LO BO ce @ Sev, SE 1 Goi €8 
GG Sy 8° TOE 6v°O ste @) ac SO°O veo yoo 8T°O 6v°O cO°V 86°0O Seat c8 
8L2 SOs 9° SOV Ee) EG AO 0°8 80°0 ve “O vo°o SGC YOO poms 8L°0 OORT 18 
SKE Ly c°TCD Som 8S°0 Gus LOO OVO vO°O CamO IEE AO) OS*? SHE 10. SSE €L 
96E Le SUES! Sia OR 90°0O ODL EO So eO ws “CO Ste NB, 00°9 66°0 S8°c CL 
PLY. S10) DP WEO) S10) = SOBOR= OGire Cleo oS) “@, 900 cee0 97 °O Gis 7S, eS SOxS BE 
BCE SOs 7 WEES) Chie 0). en On OvV-L SO0*O 8c°O €0°O vv °O Diao OS GilgaaL Oia OL 


£OOH 


fooeD)=—ss T /bu T/bu 23 TWoOTOd 23TOTeD _*os 
T/bu se uoarat SPTTOS fom’ qetl 
ssoupzeyH peaTosstaq qoedsez yqQTM uot SAZTT/SETOU-TITTW UT UOTJeAQUEeDUOD TeIOL 


TeqaoL -einjes jo verzbeq 


= = + + = ote ++ 


.s : 
? i: sf = . y : , 
‘a oO ve - 
Coe eee 


es 


WATER RESOURCES REPORT 10 


Maps and Figures in Pocket: 


Map 1—Topography, Drainage and Hydrometric Stations 

Map 2—Surficial Geology 

Map 3—Locations of Wells 

Map 4—Bedrock Topography 

Map 5—Overburden Thickness 

Map 6—Soils 

Map 7—Piezometric Surface for Bedrock 

Figure 2—Generalized Bedrock Geology, Southwestern Ontario 

Figure 3—Geologic Cross-Section of Bedrock and Overburden Formations 

Figure 5—Geologic and Hydrogeologic Cross-Sections 

Figure 14—Streamflows in the Blue Springs Creek Basin 

Figure 15—Locations of Streamflow Stations, Observation Wells, 
Precipitation Stations, Piezometers and Sub-Basins 

Figure 44—Chloride and Sulphate to Bicarbonate Milli-Equivalents 
per Litre Ratios in Water Samples 

Figure 45—Degree of Saturation with Respect to Dolomite and Calcite 

Figure 46—Calcium to Magnesium Ratio of Water Samples 


Note: Figures 14/15 and figures 44/45/46 are printed together on two sheets 
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Hydrometric station, monthly readings with station number 


Topographic contour, interval 25 feet 
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Cartography by R. Zimmermann and C. Lochan 


Base map derived from 1:50,000 sheets of the National Topographic series. 
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Pee | Eramosa member: bituminous dolomite 
Re Wiarton member: bright grey mottled dolomite 


Note: The Colpoy Bay member probably outcrops along the middle portion 
of Blue Springs Creek, but was not identified in the basin. 


SYMBOLS 


p—~ Geological boundary, approximate 
[ s\ | Measured bedrock strike and dip in degrees 


= Line of section shown in figure 3 


SOURCES OF INFORMATION 


Surficial geology compiled by J. Coward from the following references: 
Karrow, P.F., 1968, Pleistocene geology of the Guelph area; Ontario 
Department of Mines, Map 2153. 


Sanford, B.V., 1969, Geology of the Toronto-Windsor area; Geological 
Survey of Canada, Map 1263A. 


Cartography by R. Zimmermann and C. Lochan. 
Base map derived from 1:50,000 sheets of the National Topographic series. 
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Figure 14. Streamflows in the Blue Springs Creek basin used to determine the losing and gaining 
reaches along the creek. 


Government 


LEGEND 


Basin and sub-basins 


Blue Springs Creek basin 
Federal sub-basin 

Boy Scouts camp sub-basin 
Eden Mills sub-basin 
Knatchbull sub-basin 


Cedar Grove sub-basin 


Streamflow gauging station 
Precipitation station 


Observation well 


Piezometer 


Scale 
(@) 


0 


Figure 15. Locations of streamflow stations, observation wells, precipitation stations, piezometers 
and sub-basins in the Blue Springs Creek basin. 
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SOURCES OF INFORMATION 


Bedrock topography by J. Coward, 1976, on the basis of water well 
records on file with the Ontario Ministry of the Environment. 
References: 

Karrow, P. 968, Pleistocene geology of the Guelph area; Ontario 
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Department of Mines, Map 2153, 

Sanford, B.V., 1969, Geology of the Toronto-Windsor area; 
Geological Survey of Canada. Map 1263A Gore 
Cartography by M. Lakin and C. Lochan. Publications 


Base map derived from 1:50,000 maps of the National Topographic 
series. 
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Soil boundary, approximate 
Burford and Fox loams: developed from outwash sands and gravels; 
well drained 


Colwood loam: developed trom silty sand of lacustrine origin; 
c poorly drained 


Dumfries loam: developed from coarse textured till material usually on 
moraines; well drained and often containing many stones and boulders 
Don 


inybrook loam: developed from kame and esker sands and gravels; 
well to excessively drained 
Farmington loam: developed from very thin till over bedrock 
(e, Guelph loam: devel; from tills in till plains or gently sloping 
moraines; well drait and @ good crop soil 
Killean loam: imperfectly drained members of the Dumfries loam 
Lily loam: poorly drained members of the Dumfries loam 


Marsh, mesisol and peat: poorly drained swamp deposits, developed over a 
variety of parent materials; often indicates ground water discharge areas 


Parkhill and London loams: imperfectly and poorly drained members of the 
Guelph loam 


SOURCES OF INFORMATION 


Soil configuration compiled by J. Coward from the following references: 


fellington County; 
ulture, Report 35. 
ntario 
43 


Hoffman, D. W., Matthews, B. C. and Wicklund, R, E., 1963, Soil survey of We 
Canada Department of Agriculture and the Ontario Department of Agricul 
Gillespie, J. E., Wicklund, R. E.and Miller, M.H., 1971, Soils of Halton County; O1 
of Agriculture and Food and the Canada Department of Agriculture, Report 


Department 


Cartography by M, Lakin and C. Lochan 
Base map derived from 1:50,000 sheets of the National Topographic series. 
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SOURCES OF INFORMATION 


Piezometric surface by J. Cowerd 1976 on the basis of water well 
records on file with the Ontario Ministry of the Environment, 
Cartography by M. Lakin and C. Lochan 


Base map derived from 1:50,000 sheets of the National Topographic series. 
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Figure 2. Generalized bedrock geology, southwestern Ontario, after Hewitt (1964, 1972). 
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The Amabel group was tentatively split into 


three be 
well log information on the dolomite colour, and on the 
of each member: 


Eramosa member: brown bituminous dolomite 


Wiarton member: light grey, grey mottled dolomite 
Colpoy Bay member: white to blue-grey dolomite 
Reynales formation: brown dolomites and shales 
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SYMBOLS 
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Water well number, location on Map 3 


Oil and gas well number, location on Map 3 


SOURCES OF INFORMATION 


Geology by J. Coward, 1976, on the basis of water well records on file with the Ontario 
Ministry of the Environment, and oil and gas well logs assembled by the Ontario Ministry 
of Natural Resources. 


References: 
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Figure 44. Chloride and sulphate to bicarbonate milli-equivalents per litre ratios in water samples. 
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Figure 45, Degree of saturation with respect to dolomite and calcite. 
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Figure 46. Calcium to magnesium ratios of water samples. 
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